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I. INTRODUCTION 


JADEITE is a mineral to which much interest attaches by reason of the use 
which has been made of it since ancient times—especially in China—for 
the fabrication of works of art which have been greatly admired. It is a 
much rarer material than nephrite which is another mineral product to which 
the name jade is also commonly given; happily, the two substances are 
readily distinguishable from each other. In its finest forms, jadeite is known 
as gem-jade or imperial jade and differs from the commoner varieties in its 
optical characters; it exhibits a beautiful green colour resembling that of 
emerald and is semi-transparent instead of being cloudy or opaque. There 
is an extensive literature which deals with jadeite from the mineralogical 
point of view; works on gemmology also naturally devote a good deal of 
space to it. We have not however found in the literature any explanation 
of the optical characters which jadeite exhibits in terms of its structure. 
Indeed, it does not appear to have been realised that there is a problem 
awaiting answer in that connection. The opportunity for examining the 
subject presented itself to us when a collection of specimens of Burmese jade 
which was acquired for the Museum.of this Institute became available for 
our studies. 


2. THE STRUCTURE OF JADEITE 


The chemical composition of jadeite corresponds to the formula 
NaA\(SiO;), and indicates that it belongs to the series of monoclinic pyro- 
xenes of which diopside CaMg(SiO;), is a typical example. X-ray studies! » % 
support the latter conclusion, the X-ray powder diagrams of jadeite and 
diopside showing a close resemblance to each other. The published 
chemical analyses‘ show the presence in jadeite of varying proportions of 
other oxides, notably the oxides of calcium and magnesium, thereby indicat- 
ing that diopside is present as a minor constituent in the mineral. Physically, 
jadeite is a polycrystalline aggregate. Examination of thin sections under 
the polarising microscope reveals the presence of individual crystals which 
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are birefringent, their greatest and the least refractive indices being 1-667 
and 1-654 respectively. 


Whon light enters a polycrystalline aggregate in which the individual 
grains are birefringent, it would necessarily suffer reflection and refraction 
at the intercrystalline boundaries. The smaller the individual crystallites, 
the more numerous would be such reflections and refractions. Hence, on 
this view, none of the incident light would succeed in penetrating through 
the material, if it is fine-grained and also of sufficient thickness. In the 
absence of any absorption properly so-called, the incident light would mostly 
be sent back towards the source by repeated internal reflections. In other 
words, a fine-grained aggregate of birefringent crystals should present a 
high albedo by reflected light and appear opaque in transmission. 


The foregoing inferences based on geometrical optics do not however 
represent what is actually observed in the case of jadeite. Indeed, the con- 
verse is much nearer the truth; in other words, the specimens exhibiting 
a coarse structure are those which are opaque to light and reflect it strongly, 
while the specimens which are most readily penetrated by light are those 


having a fine-grained structure. Why this is so is the problem considered 
in the present paper. 


3. THe X-RAY DIFFRACTION PATTERNS 


To cstablish the relationship between structure and optical characters 
briefly indicated above. we have recorded the X-ray diffraction patterns of 
specimens exhibiting diverse optical behaviour. Three of these diagrams 
which are representative of the rest are reproduced in Plate XII. The film 


to specimen distance wa; in every case 5 cm. and MOKa radiation was 
employed to record the patterns. 


Figure | was recorded with a piece of jade about half a millimetre thick 
cut out from one of our specimens and polished flat on both sides. The 
X-ray pencil passed normally through the specimen. The plate exhibited 
a notable measure of translucency; when placed on a printed sheet of paper, 
the finest print could easily be read through it. The X-ray record shows 
clearly that the crystallites present in the specimen are of varying sizes and 
that crystallites of very small sizes are the most numerous; though individual 
bright spots appear in the pattern, there are also numerous fainter ones which 
run into each other and appear as continuous rings. 


Figure 2 is an X-ray diagram taken with similar arrangements, the X- 
ray pencil now grazing the edge of a flat ellipsoidal piece of green jade. It 
will be seen that the pattern is of the same general nature as that reproduced 
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in Fig. 1 except that the continuity of the rings is less evident and the 
spottiness more conspicuous. 


Figure 3 is the pattern recorded with a greenish-white specimen which 
was nearly opaque to light. It will be seen that the pattern consists of a 
number of individual spots, some of which are very bright. But hardly 
anything in the nature of continuous rings can be discerned. In other words, 
the specimen consists of relatively large crystallites. 


4. THEORETICAL CONSIDERATIONS 


To understand the optical behaviour of jadeite, we have to put aside 
the considerations based on geometrical optics and deal with the matter 
from the standpoint of the wave theory of light. It is useful, in the first 
instance, to think of an extreme case in which the individuals in the poly- 
crystalline aggregate are assumed to be of very small dimensions even in 
comparison with the wave-length of light and are also randomly orientated. 
Such an aggregate could obviously be regarded as equivalent to an optically 
isotropic medium with a refractive index which is the average of the principal 
refractive indices of the crystallites. A beam of light incident normally 
on a slab of such material would traverse it as a coherent wave-train, and 
emerge from the rear surface of the slab, but not the whole of the incident 
energy would succeed in so emerging. For, in consequence of the varying 
orientations of the crystallites and their birefringence, the medium would 
exhibit local variations in optical behaviour, and there would therefore be a 
diffusion of light both forwards and backwards. In the limiting case when 
the grains are very small, the diffusion would extend approximately with 
equal strength both forwards and backwards. More generally, however, 
there would be a concentration of the diffused radiation in the forward 
directions surrounding the path of the regularly transmitted light-b:am. 
There would also be a diffusion of light backwards towards the light-source, 
but this would be less conspicuous. 


The situation set forth above would be modified if the grains of the 
medium are of dimensions comparable with the wave-length of light. Never- 
theless, provided the birefringence of the individual crystallites is small, 
we may continue to regard the aggregate as an isotropic medium which 
exhibits local variations in refractive index depending on the direction of 
vibration in the light beam traversing the material. These variations would 
result in the phase of light-waves passing through an individual grain being 
either accelerated or retarded as the casc may be. The errors in phase thus 
resulting when summed up over the total path would tend to cancel each 
other out; the residual errors in phase would be the smaller relatively to 
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the total path, the more numerous the particles are, in other words the 
smaller is their size. Thus, it follows that if the birefringence is not large 
and the material is very fine-grained, the waves of light would succeed in 
penetrating the medium. Nevertheless, the residual errors in phase would 
result in a reduction of the absolute intensity of the transmitted beam and 
ultimately also its total extinction. The missing energy would appear for 


the most part in a cone of diffuse radiation surrounding the direction of 
propagation of the incident light-beam. 


Wave-theoretical considerations also indicate that the proportion of 
the incident energy sent back towards the source as a diffuse reflection by a 
fine-grained aggregate would be far less than that indicated by considerations 
based on geometrical optics. For, the phase of a reflected wave relatively 
to the incident wave depends on whether the first or the second medium on 
the two sides of the reflecting boundary has the higher index. The bounda- 
ries encountered by a light-wave traversing a polycrystalline aggregate would 
obviously be of both kinds and in approximately equal numbers. The 
reflection in the two cases being in opposite phases, their effects would tend 
to cancel each other more or less completely in all directions. 


Thus we may sum up the situation by the statement that light would 
penetrate into a polycrystalline aggregate much more freely and would be 
reflected backwards much less freely than that might be anticipated in terms 
of geometrical optics, provided that the aggregate is sufficiently fine-grained. 


Thus, we have an explanation in general terms of the semi-transparency 
exhibited by the finest varieties of jadeite. 


5. SUMMARY 


The phenomenon of semi-transparency exhibited by the finer varieties 
of jadeite is discussed. X-ray diffraction studies and theoretical considera- 
tions based on wave-optics alike indicate that this is a consequence of the 
very small particle size of the crystallites forming the aggregate. 
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1. 


CapMIUM sulphate (3CdSO,.8H,O) differs from most other crystalline 
sulphates containing water of crystallisation in that the proportion of water 

is very much less and the crystal structure analysis (Lipson, 1936) shows 

that the forces between the cadmium and SO, ions are so large that to a first 
approximation one can regard the structure as composed of Cd and SO, 
groups in which the water molecules occupy the interstices. Although the 
spectrum of its water of crystallisation has been studied and compared with 

those of other crystals containing water of crystallisation (Hibben, 1937; 

P. G. N. Nayar, 1938) its Raman spectrum has not been so far investigated 

in any great detail. Its Raman spectrum recorded by Hibben and reproduced 

in his book ‘ Raman Effect and its Chemical Applications ” (1939) does not 

show much detail and Hibben limited himself to pointing out the splitting 

in the crystalline state of the triply degenerate oscillation (1106cm.) of 

the free SO, ion and no further study of the spectrum was made by him. 
Besides these, the strongest Raman line 1003 cm.-? alone was recorded by 
Embirikos (1930) and compared with those of other sulphates. An investiga- i 
tion of its Raman spectrum made by us recently revealed such an abundance 

of detail due to the use of spectrographs of large dispersion, improved Rasetti 
technique and more efficient filtering of the exciting radiation that it was 

felt desirable to report the results given below. 


INTRODUCTION 


2. EXPERIMENTAL DETAILS 


The crystals used were prepared from the aqueous solution of the salt 
by the method of slow evaporation and were about 1 cm.x5 mm.x5 mm. 
in size with well-developed faces. Besides those grown in our laboratory, 
one other fine specimen also grown from solution by the same method was 
made available to us through the !'‘ndness of Prof. J. P. Mathieu to whom 


our sincere thanks are due, 
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The resonance radiation of mercury (A 2537) from a specially designed 
water-cooled quartz arc placed in the field of an electromagnet was used to 
excite the Raman spectrum. Using the medium quartz spectrograph and 
slit width of -025 mm. exposures of the order of 2 hours were found to be 
sufficient to record the complete spectrum. It is worth mentioning here 
that though the experimental arrangement was found to be adequate to 
record the symmetric oscillation (1006 cm.-') of the SO, ion and the water- 
bands in a couple of minutes, even exposures of the order of 2 days failed 
to bring out the second order lines. In order to study the lattice spectrum 
in more detail a large quartz spectrograph with a dispersion of about 40 cm.-! 
per mm. in the A 2537 region was employed. With this instrument using 
a slit width of -05 mm. an exposure of the order of a day was enough to 
record the complete Raman spectrum. 


3. RESULTS AND DISCUSSION 


The Raman spectrum of cadmium sulphate 3CdSO,.8H,O consists of 
23 frequency shifts all of which except the waterbands are marked in Fig. 1 b 
of Plate XIII. Fig. 1 (a) is the spectrum taken with the large quartz spectro- 
graph and enlarged to the same extent as its microphotometer record shown 
in Fig. 1(b). The waterbands as well as the other Raman lines could also 
be easily identified in Fig. 2 a which is an enlargement of the spectrogram 
taken with the smaller instrument while Fig. 2 5 shows the direct mercury 
spectrum. The observed frequency shifts are 45 (2), 56(7), 64(7), 75 (2), 
89 (2), 104 (5), 124 (4), 138 (5), 170 (2), 190 (3), 219 (5), 261 (1), 331 (3), 
456 (6), 462 (5), 626 (2), 1006 (30), 1077 (3), 1101 (9), 1115 (10), 1157 (3), 
3347 (7) and 3428 (15) cm.-! the figures within brackets representing the 
visually estimated relative intensities. The values of the frequency shifts of 
all the lines except those of the waterbands were determined from measure- 
ments made on the spectrogram taken with the El (large) spectrograph and 
may therefore be regarded as accurate. We will now consider the observed 
spectrum in relation to the known crystal structure. 


The crystal structure of 3CdSO,.8H,O has been determined by Lipson 
(1936) and reveals certain interesting features. Unlike most other sulphates 
with water of crystallisation, like alum, KAI(SO,),12.H,O or Epsom salt 
MgSO,-7H,O where the water is known to form an essential part of the 
structure, in cadmium sulphate eight-third hydrate the main structure is 
governed by the forces between the metal atoms and the SO, groups and 
the water molecules just fill up the interstices. The crystal belongs to the 
monoclinic space group I2/h (C*%,) and contains four molecules of 
3CdSO,.8H,0O in the unit cell. In the crystal, only 4 Cd and 4 SO, occupy 
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special positions and are located on the two-fold axis and all the remaining 
atoms are in general positions. The local symmetry of even the special SO, 
group in the crystal is only C; and as such it is only reasonable to suppose 
that on account of this and also the known multiplicity of ions in the unit 
cell, the degeneracies of the frequencies belonging to tetrahedral SO, ion 
areremoved. Actually one observes only 8 frequencies which could be attri- 
buted to the internal oscillations of the SO, ion as shown below. 


RAMAN FREQUENCIES OF THE SO, ION IN cm. 


In the free state> 454 622-983 1106 
| | | 


In 3CdSO,.8H,O- 456 462 626 1006 1076, 1101, 1115, 1157 


The Raman line 1076 cm.-! has been explained by Hibben as a combination 
of the frequencies 450 and 626 cm.!; this appears to be erroneous since in 
the crystal the line 1076 cm.“ is very intense while the fundamental 626 cm.“! 
itself is very weak. While the degeneracy of the oscillation (1106 cm.-') 
of the free SO, ion is seen to be removed completely and splitting is evident, 
one notices, however, a different behaviour of the other two degenerate 
oscillations. In the crystal the Raman lines due to the splitting of the doubly 
degenerate oscillation are neither well resolved nor sharp as for example 
in gypsum which also belongs to the monoclinic class and is known 
to contain loosely-bound water of crystallisation (Wooster, 1936). Also in the 
present study only a very weak line corresponding to the triply degenerate 
oscillation of the free ion is observed at 626cm.- while in gypsum three 
fairly sharp lines (618, 622 and 672cm.—') have been recorded by R. S. 
Krishnan (1945) in the same region. Further the spectrum of the water of 
crystallisation is comparatively simple and consists of only two bands the 
frequency shifts of which agree well with those reported by Hibben (1937) 
and by P. G. N. Nayar (1938). These observed features of the lines due 
to the internal oscillations of the SO, ions and of the water of crystallisation 
suggest that in the crystalline state while the influence of the additional forces 
due to the neighbouring ions is by no means negligible, it does not lead to 
any major distortion of the disposition of the bonds from their idealized 
nature as expected by Lipson (loc. cit.). 


Besides the 8 Raman lines due to the internal oscillations of the SO, 
groups and the two waterbands, the Raman spectrum contains 13 low fre- 
quency shift lines which arise from lattice oscillations. Neglecting the water 
molecules, a simple group theoretical calculation shows that the number of 
translatory type of oscillations permitted in Raman effect is 36, of which 
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16 are of type A,g and the remaining 20 of type B,g. However, we find that 
even under the high dispersion of the large quartz spectrograph the crystal 
exhibits only 13 Raman lines which are attributable to the translatory type 
of lattice oscillations. Since it is known that for this crystal only the 
vibrations of type Ayg will show a variation of intensity with orientation an 
attempt was made to study the variation of the intensity with orientation 
with incident unpolarised light. The observations indicate that the lines 
219, 104, 138 cm.-! may tentatively be assigned to type A,g and those with 
frequency shifts 170, 190, 331 cm.-! to the type Bg. It was also observed 
that among the Raman lines due to the SO, ions and the waterbands those 
at 462 and 3347 are of type Axg. 


It is a pleasure for us to thank Professor R. S. Krishnan, for his help 
and encouragement. 


SUMMARY 


The Raman spectrum of a single crystal of cadmium sulphate 
(3CdSO,.8H,O) has been investigated using the resonance radiation of 
mercury as exciter. The recorded spectrum reveals 23 Raman lines and of 
these 8 have been attributed to the internal oscillations of the SO, ion and 
2 to the oscillations of the water molecules. While the number of trans- 
latory type of lattice oscillations permitted by the selection rules is 36, only 
13 are observed experimentally. A tentative assignment for some of these 
vibrations has been made from a study of the variation of intensity of the 
lines with orientation. 
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FIG. 2. (a2) Raman spectrum of crystalline cadmium sulphate taken with the medium 
quartz spectrograph. 


(4) Direct Mercury Spectrum. 
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1. INTRODUCTION 


SouND absorption in most liquids is proportional to the square of frequency 
as required by the classical formule for viscosity and heat conduction losses. 
Departures from this law are now generally ascribed to relaxation phenomena 
in intermolecular processes. Herzfeld and Rice (1928) followed by Kneser 
(1948), Markham (1950), Liebermann (1949), Hall (1949) and Lamb and 
Pinkerton (1949) have offered explanations based on relaxation. Herzfeld 
and Rice (Joc. cit.) attributed anomalous absorption to what may be de- 
scribed as thermal relaxation, that is a lag in the adjustment of energy 
between the external and internal degrees of freedom of the individual mole- 
cules during passage of a sound wave. Kneser, Markham and others have 
treated the problem more generally and visualised other types of relaxation. 
In associated liquids there may be changes in the degree of association in 
the presence of a sound wave, which may take a finite time and hence lead 
to absorption of energy. Liebermann postulated a disturbance in the 
dynamic equilibrium between the solute and its ions in an aqueous solution 
of MgSO,, to explain the observed peak in absorption at a particular 
frequency in such solutions. Whatever may be the relaxation mechanism 
proposed, some common conclusions follow from all these theories. One 
is that there is a frequency w,, called the relaxation frequency, and an 
associated relaxation time 7 characteristic of each liquid. In regions not 
too close to this frequency, absorption is linearly related to the square, and 
velocity is independent, of frequency. On the other hand, near this fre- 
quency, velocity shows dispersion. At a frequency which is a simple func- 
tion of w, and less than the latter, the absorption per wave-length (aA) goes 
through a peak when plotted against frequency. Sulphur dioxide, ammonia 
and carbon dioxide have all been studied thoroughly in their gaseous states 
and are known to exhibit the abovementioned relaxation phenomena. It is 
considered interesting to see how they behave in solution. In the present 
work, solutions of these gases in some typical liquids are studied for the 
first time, 
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2. EXPERIMENTAL METHOD 


Employing a method similar to that described by Pellam and Galt 
‘ (1946), in which instead of a reflector, a second crystal is used as the receiver, 
sound absorption measurements have been made in several gaseous solutions. 
The two crystals used have identical fundamental frequencies of 2-12 Mc/s 
and are worked at their 3rd, 5th and 7th harmonics besides their funda- 
mental. The temperature is controlled to within + 0-5°C. The alignment 
of the crystals as well as the accuracy. of the method as a whole is tested by 
taking measurements in double distilled and suspension—free water at all 
the frequencies. For attenuation measurements, the height of the received 
pulse is adjusted to be constant by using a calibrated standard decade 
attenuator with a bandwidth of 0 to 30 Mc/s and capable of reading correct 
to 0-1 db. The attenuation measurements are accurate to within 3 to 6 per 
cent., the lower limit applying to the higher frequencies. 


Sulphur dioxide used in this work is manufactured by I.C.I. as a 
refrigerant and is of guaranteed purity. Ammonia is obtained as an aqueous 
solution (sp. gr. 0-88) of analytical purity. Carbon dioxide is a commercial 
sample obtained locally. Sulphur dioxide solutions are estimated by titra- 
tion in the usual way. Strengths of ammonia solutions are determined 
from density data. Only a saturated solution of carbon dioxide has been 
studied and no attempt is made to determine its strength. 


3. RESULTS AND DISCUSSION 


The highest concentrations of sulphur dioxide studied are 0-1M in 
water, 0-55 M in Benzene, 0-95 M in Ethyl Alcohol and 0-3 M in Nitro- 
benzene. Ammonia is dissolved only in water and the maximum concentra- 
tion studied is 0-55M. All these solutions are investigated at the room 
temperature. Carbon dioxide which is very sparingly soluble is examined 
at the low temperature of 8° C. so that more of the gas will dissolve in water. 


A representative set of results for all these solutions are presented in 
Figs. 1 and 2. 


From the above graphs it is seen that solutions of carbon dioxide and 
ammonia in water behave normally in that a A increases linearly with frequency 
at all concentrations. Therefore solutions of these gases in other solvents 
are not studied. On the other hand, aqueous sulphur dioxide solutions 
show a peculiar behaviour, in several ways. The attenuation increases very 
rapidly with concentration. The plot of ad against frequency is not a 
straight line with a positive slope. It actually exhibits peaks which appear 
to occur at lower frequencies for lower concentrations, At the higher 
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concentrations, the frequency of peak absorption (aA) reaches a limit and 
becomes independent of concentration. All these peculiarities point to a 
relaxation phenomenon and the characteristic relaxation frequency at the 
higher concentrations lies in the range 2 to 15 Mc/s while at lower con- 
centrations it is less than 2 Mc/s. As a probable cause for this phenomenon, 
one can think of the gas itself or the chemical reaction going on in the 
aqueous solutions. That the gas is not the cause of this phenomenon 
follows from the fact that non-aqueous solutions of sulphur dioxide do not 
show any such relaxation. This leaves the dynamic chemical equilibrium 
between H,SO, and H+ and HSO; ions, namely H,SO; = Ht + HSO,, 


¢ SO, in Ethyl! Alcohol 
SO. in Benzene 

1504 * SO. in Nitrobenzene 
@~----e-~--~-e NH; in water 
CQ, in water 
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as the probable cause. In the presence of the ultrasonic beam, especially 
n the particular range of frequencies used in this work, there seems to be a 
change in this equilibrium followed by absorption of energy. This means 
that the above reaction is proceeding in either direction at a rate comparable 
with the period of oscillation of the sound waves and that with increasing 
concentration it is increasing upto a limit. 


Studies in aqueous solutions of bisulphites of potassium and sodium 
as well as sodium sulphite show similar results. This confirms that where 
‘ an HSO, ion is concerned, relaxation phenomena are again exhibited. 
Details of these investigations, which are of a confirmatory type, will be 
published separately. 


It is, however, necessary to explain why aqueous solutions of ammonia 
and carbon dioxide which also undergo chemical changes in solution, do not 
exhibit relaxation effects. These three gases have been thoroughly studied 
by Fricke (1939) and Petralia (1953) among others and they have found 
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that carbon dioxide exhibits a peak in absorption at about 20Kc/s. The 
corresponding values for sulphur dioxide and ammonia are 1-040 Mc/s. 
and 393 Mc/s. respectively. The relaxation times in the latter have been 
calculated as 1-85x10-7 sec. and 8-1x10-™sec. respectively. It means 
that these are the times in which the extra translational energy supplied by 
the sound beam is converted by the respective molecules into internal energy. 
It is reasonable to expect that the conversion of chemical energy into kinetic 
energy or back also involves times of this order of magnitude for the two types 
of molecules. Thus, the frequency range to which the present studies are 
confined explains the occurrence of a relaxation phenomenon in the case 
of sulphur dioxide solutions only and not in ammonia solutions. On the 
same basis, the rate of reaction in the case of carbon dioxide will be so slow 
that so far as the frequencies involved in the present experiments are con- 
cerned, it is as good as not taking place at all. 


4. SUMMARY 


Absorption of ultrasonic energy by solutions of sulphur dioxide ammonia 
and carbon dioxide in water, and in the case of the first gas in several other 
organic solvents, has been studied by the pulse technique. Aqueous sulphur 
dioxide solutions exhibit a relaxation phenomenon in the frequency range 
of 2 to 15 Mc/s. as evidenced by the absorption (a A) showing a peak in this 
region when plotted against frequency. This peak occurs at higher frequencies 
for higher concentrations upto a limit after which the peak absorption fre- 
quency does not change with concentration. These observations, along 
with the fact that no such effect is observed in the case of other solvents, 
lead to the conclusion that the phenomenon is caused by a change in the 
dynamic equilibrium between H,SO,; molecules and H*+ and HSO;- ions 
in the presence of the sound beam. The absence of a similar effect in the 
case of ammonia and carbon dioxide solutions which also involve a chemical 


reaction is explained on the basis of the widely differing relaxation times of 
the three gases. 


In conclusion, the author acknowledges with grateful thanks the guidance 
and help given by Professor S. Bhagavantam who suggested this problem. 
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1. INTRODUCTION 


WHEN polarised light of any form is incident on a birefringent plate, the 
phase retardation 6 introduced between the two waves during their passage 
through the plate, is far from being the same for all wave-lengths; in fact, 
being given by 5 = (27/A) (u, — 12) ¢, this retardation of phase varies almost 
inversely as the wave-length A, that is, if the birefringence (u, —p,) does not 
disperse notably with wave-length. So much so, that what acts as a quarter- 
wave retardation plate for the deep red end of the spectrum will, for the 
wave-lengths in the deep violet, behave practically as a half-wave plate. 


The question of the achromatisation of devices in this field has engaged 
the attention of several workers. References to their investigations’ may be 
found cited in the paper entitled, “‘ Réalisation d’un quart d’onde quasi-achro- 
matique par juxtaposition de deux lames crystallines de méme nature” by 
Destriau and Prouteau!; the two plates referred to, do not, of course, have 
their principal planes parallel, and the ‘compound plate’, comprising the 
two superposed plates, can transform incident circularly polarised light to 
plane polarised light vibrating at a certain azimuth, or vice-versa. But the 
device described by these authors would more properly be called an achro- 
matic circular polarizer (or analyser); to:call it an achromatic quarter-wave 
plate would be incorrect since the combination does not have the usual 
at ributes of a quarter-wave plate—it cannot, for example, be used for the 
analysis of elliptically polarised light in the usual manner of an ordinary 
retardation plate. 


Nevertheless it is as a circular polarizer or analyser that a quarter-wave 
plate is often used—as when it is inserted with its principal planes at an angle 
of 45° to those of a nicol, in a petrographic microscope. And we shall 
reserve for the second part, the problem of superposing birefringent plates 
in such a manner that the combination as a whole behaves as an achromatic 
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quarter-wave plate. In this paper we shall describe a circular polarizer 
which is expected to have a much higher degree of achromatism than the one 
discussed in the paper quoted above. The description of the device may be 
found at the end of the next section; its applicability is not limited to the 
visible. The birefringent plates in our discussion are all of the same material ; 
the dispersion of the birefringence of this material is not assumed to be 
negligible—though such an assumption would indeed be justified for plates 
of muscovite mica, at least for the visible and ultra-violet wave-lengths.* 


2. THE ACHROMATIC CIRCULAR POLARIZER 


The Poincaré sphere* lends itself very conveniently for the theoretical 
discussion of the two problems at hand—especially the one to be treated 
later. (As is well known, any point on the surface of this sphere, of latitude 
2w and longitude 2/, represents an elliptic vibration of ellipticity | tan w| 
the major axis of which makes an angle / with a fixed reference direction. 
Passage through a birefringent plate of retardation 5, the orientation of 
whose principal planes are given by O and O’, is equivalent to the operation 
of rotating the sphere by an angle 5 about th> equatorial diameter OO’.) 


=} 


A quarter-wave plate with its fast vibration direction in the orientation 
O; (Fig. 1) can transform the incident linear vibration R—inclined at 45° 
to the principal planes—to circularly polarised light represented by the 
pole S—but only for the wave-length A for which its retardation is exactly 


* Expositions in English, of the properties of the Poincaré sphere are apparently not so 
common in text-books, but may be found in numerous recent articles: see e.g., Ramachandran 
and Ramaseshan, J. Opt. Soc. Am., 1952, 42, 49. 
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m/2. For the wave-lengths A, and A, on either side of A, for which the respec- 
tive retardations are 7/2 — « and a/2 + e, the emergent light will be cir- 
cularly polarised only if the incident vibrations are given respectively by 
R, and R,; these two points lie on the meridian SR and are equidistant 
from R (arc R,R= «). Thus if the ellipticity of the incident vibrations 
were dispersed roughly along R,R. in a proper fashion, the quarter-wave 
plate would transform the vibrations for all the wave-lengths to -approxi- 
mately circular polarised light. 


The required dispersion of the ellipticity of the vibrations incident on 
the quarter-wave plate, it may be seen, can be effected by allowing plane 
polarised light P to first pass successively through two half-wave plates, 
the orientation of whose fast axes are given by O, and O, respectively and 
are to be determined. In order that the light emerging from the third 
retardation plate (i.e., the quarter-wave plate) should be exactly circularly 
polarised for all the three wave-lengths A,, A and A,, the state of polarization 
of the light incident on the second half-wave plate for these three wave-lengths 
should be given by Q;, Q and Q, respectively—got by the following 
construction. The point Q is marked off on the equator such that O,.Q = O.R 
(= a say); the arc O,Q, is drawn such that R,O,Q, = 7 — 2e, and O,Q, 
= O.R, (= b say); Q, is a point symmetrically placed with respect to Q, 
on the lower hemisphere. The first half-wave plate has to be oriented in a 
position depending on the quantity a which determines the orientation of 
the second; for we must have 0,Q, = 0,Q, = 0,Q (= e say); and further 
we must choose a such that PO,Q, = 7 — 2e. If the plane of the incident 
vibration P be now set such that O,P = e, then the light emerging from 
the combination will be circularly polarized for the three wave-lengths A, 
A and A, for which the retardations of a quarter-wave plate would be 7/2 —e, 
m/2, and 7/2 + « respectively; ¢ is an arbitrarily chosen parameter on which 
will depend the two opposing characteristics of the combination: the range 
and the degree of its achromatism. 


By the use of spherical trigonometry we shall get an equation for a 
involving the parameter «. Thus, from the equilateral triangl> Q,0,Q, 
denoting arc Q,Q by c 


cos 2e = — cos? Q,QO, + sin? Q,QO, cos c 
1 + cos 2e = 2sin? Q,QO, cos? c/2 
From the triangle Q,0.Q we have, writing Q,O.Q = C, 


or 


sinc _ sind 
snC sin Q,QO, 


Achromatic Combinations of Birefringent Plates—I 


i.e., 
c sin b sin C 


We thus obtain 
sin? b sin? C = 2 sin* 5 ( + cos 2e) (1) 
We shall express the quantities on the left-hand side in terms of a. Denoting P 
R,0,.R by C’ 
sin C = sin (2e — C’) = sin 2e cos C’ — cos 2¢ sin C’ (2) 


From the right-angled triangle R,O.R 


and 


osb i 
cos = cos a cos - == COS sin @ 
sin a sin b 


since 
cos b = cos e cosa 


Substituting the expressions obtained for cos C’ and sin C’ in (2) 


sin € 
sin C = ——, (2 cos? e sin a — cos 2e) 
sin b 


Introducing this in (1) 
(2 cos? « sin a — cos 2)? = 2 cot? «(1 — cosc) (4) 


It remains to express cosc in terms of a. From triangle Q,QO, 


cos c= cosacosb+ sina sinbcosC 
cos C = cos (2« — C’) = cos 2e cos C’ + sin sin C’ 


cose 
sin 5 


Usin_ this relation, together with (3), the expression for cos ¢ becomes, 


(cos 2esina + 2sin%«) | 


cos c = cos «(1 + 2 sin? « sin a — 2 sin? ¢ sin? a) 


Sutstituting for cos c in (4), and rearranging, we finally get an equa‘ion in 
sina of the form: 


Asin?ta+ Bsinga+C=0O (I) 


A2 


sn C’ = si b 
‘ 
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A = 4cos® e (1 — cos e) 


B = 4cos? « (cos 2e — cos e) 

_ C= 2cot* e(1 — cos e) — cos? 
An expression for e may be obtained from triangle Q,0,Q 
cos c = cos* e + sin* e cos 


~ —cos2e 
Using (4), 
2 cos* sin a — cos 


sine = (Il) 


The combination can produce exactly circularly polarized light for three 
wave-lengths. If we neglect the dispersion of the birefringence, the two 
extreme wave-lengths will obviously be (1 + /) A and (1 — f) A—where we 
have denoted the ratio «:7/2 by f; the achromatism of the combination 
may however be considered to extend up to two wave-lengths outside this 
interval, for which the deviation from circular polarisation is roughly the 
same as the maximum deviation attained within this interval—which latter 
may be expected to occur in the neighbourhood of the wave-lengths 
(1 + -5f)A. Thus we may with some arbitrariness, consider the combina- 
tion to be achromatic within the range (1 + 1-25/)A. 


Choosing 2¢« = 47° (which means 1-25 f= -325), equation (I) when 
solved gives sina = +6642 or a= 41° 37’. 


The corresponding value of e obtained from (II) is e = 13° 43’. 


On the other hand, if we choose 2« = 36° (which means 1-25 f= -25) 
we obtain 
a = 43° 8’; e = 13°. 


Thus to summarize our results, an achromatic arrangement for producing 
circularly polarised light is obtained by allowing parallel plane polarised 
light to pass normally through a combination of three superposed plates of 
the same material; the first two should be half-wave plates, and the last, 
a quarter-wave plate for the wave-length A in the centre of the spectral range 
to be covered. Let 6 be the angle by which the fast vibration direction of 
the first plate is turned with respect to the azimuth of the incident vibration; 
6,, the angle by which the fast axis of the second plate is turned with respect 
to the first; and 6, the corresponding angle between the fast axes of the 


where 
or 
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third and second plates—all the angles being measured in the same sense. 
The set of values 


0, = 4 (m/2 + a) = 65° 49’; 0, = 4 (a + e) = 27° 40’; e = 6°52’ 


can be used for covering the range from 1-325 A to -675 A and is therefore 
suited for covering the entire visible spectrum if we choose \ = 6000A. A 
greater degree of achromatism can be attained at the expense of restricting 
the range to the major portion of the visible. Thus the set of values 


6, = 66° 34’; 0, = 28° 4’; 6° 30’ 


can be used for covering the range from 1-25A to -75A. If the dispersion 
of the birefringence of the material used is not negligible, the values given 
above for the range over which the combination may be considered achromatic 
would have to be altered. 


On turning the polarizing nicol to a perpendicular position, circularly 
polarized light of the opposite sense is produced. In fact in passage through 
any succession of ‘ elliptically birefringent’ plates, two orthogonal states 
of the incident polarization will correspond to two orthogonal states of the 
emergent polarization. This is a consequence of a more general fact; the 
angular separation on the Poincaré sphere of two possible states of the inci- 
dent polarization will be equal to the angular separation of the two cor- 
responding states of the emergent polarization, since angular relationships 
temain invariant under any number of rotations. Thus if the nicol had been 
turned by 45°, the emergent wave-lengths would all have been practically 
plane polarized, but with the azimuths of the vibrations dispersed—an interest- 
ing illustration of the fact that the combination does not have the properties 
of a quarter-wave plate. 


3. EXPERIMENTAL VERIFICATION 


Two half-wave plates and a quarter-wave plate, all prepared of mica, 
were cemented together (using copal varnish) with their principal planes at 
inclinations slightly different from the first set of values given above—in order 
to cover an even wider range. The combination was laid on a mirror. With 
the aid of another mirror, white light from a point source was made to pass 
normally through a polaroid and the achromatic combination, and then 
back again—after which it reached the eye. As the polaroid is rotated a 
position is reached where the image of the source is completely extinguished ; 
this is the position where a single passage through the compound plate would 
give circularly polarized light, while the double passage gives plane polarized 
light which is crossed out by the polarizer itself. When the experiment is 
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repeated with a single quarter-wave plate the image is never completely 
extinguished, but is highly coloured. 


The two A/2 retardation plates may first be cemented together at the 
proper angle @,, using the fact that the two together will rotate the plane 
of polarization of any incident linear vibration (of wave-length 2) by an angle 
a — 26,; (this can be proved by a construction given in Part II); the quarter- 
wave plate may then be cemented and adjusted, while the cement is still wet, 
to the position of best achromatism—as determined by the test given above. 


SUMMARY 


Circularly polarised light is obtained by superposing two half-wave 
plates and one quarter-wave plate, all of the same material, such that the 
fast vibration-directions of the successive plates make specific angles 4,, @,, 0, 
with the azimuth of the linear vibration incident on the first plate. The 
required range of achromatism determines the optimum values of the angles, 
Thus, using mica retardation plates prepared for Na 5890, the range from 
4000 A to 7800 A is covered with 0, = 6°52’, 6, = 34° 32’, 0, = 100° 20’; 
while the range 4400 A to 7400 A can be covered with superior achromatism 
by taking 9, = 6°30’, 6, = 34°34’, = 101° 8’. 
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1. INTRODUCTION 


A COMBINATION of birefringent plates which can be used for producing 
circularly polarised light over a fairly large range of wave-lengths has been 
described in Part I. And it has there been pointed out that the arrangement 
cannot be properly described as an achromatic quarter-wave plate. A quarter- 
wave plate when rotated between crossed nicols will show two extinction 
positions; the achromatic combination described will show none. And to 
understand under what conditions a combination of plates will have the 
properties of a single birefringent plate—let alone a plate which is in addition 
achromatic—it is desirable to have a general method for combining the 
action of a pile of birefringent plates when their principal planes are not 
parallel. We shall follow a method differing slightly from that which may 
be found described by Pockels.! 


2. GENERAL CONSIDERATIONS 


The problem of designing an achromatic elliptic polariser is the problem 
of transforming a particular state of polarization represented by P, (on the 
Poincaré sphere) to another particular state P,’. But, of an achromatic 
retardation plate we require that it transforms every P to a corresponding 
P’ obtained by rotating the sphere by an angle 5 about an equatorial diameter. 


Passage of monochromatic polarised light through a succession of bire- 
fringent plates corresponds to a succession of rotations of the Poincaré 
sphere about a series of equatorial diameters. Any succession of rotations 
of the sphere can, however, be compounded into one single resultant rota- 
tion. Hence the action of a combination of birefringent plates corresponds 
to a rotation ¢ of the Poincaré sphere about some diameter EE’—which will 
in general be inclined to the equatorial plane; this means, of course, that 
the combination is equivalent to a single plate having elliptic birefringence— 
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E and E’ being the orthogonal elliptic vibrations propagated without change 
of form, and between which a phase retardation ¢ is incurred. 


The condition that the combination should act as a purely birefringent 
plate (for a particular wave-length A) is that the axis EE’ of the resultant rota- 
tion should be on the equator; so that, if 2w denote the lattitude of E, the 
faster ellipse, 2 = 0. And if the retardation of the plate is to have a 


required value 5, the magnitude of the resultant rotation ¢ should be equated 
to this value. 


Both the magnitude ¢ as well as the axis EE’ of the resultant rotation 
will in general alter with the wave-length, since the magnitudes (though not 
the axes) of the component rotations are a function of wave-length. Hence, 
for the combination to be considered achromatic in the immediate neighbour- 
hood of the wave-length A, some more conditions have to be imposed; when 
the wave-length is increased by 4A the values of the composite phase retarda- 


tion ¢, as well as the latitude 2w and longitude 2/ of E, should remain un- 
altered at least to the first order of approximation. 


Instead of following matrix methods, we shall use a well-known geo- 
metrical construction (illustrated in Fig. 1) for finding the resultant of any 


aN 


Fic. 1. Composition of Rotations: 2R, -+ 2R, = 2R. 


two rotations. 


If ABC be a spherical triangle described on a sphere whose 
centre is O then: a rotation about AO through twice the internal angle at A, 
followed by a rotation about BO through twice the internal angle at B is equiva- 
lent to a rotation about CO through twice the external angle at C. 


3. A BIREFRINGENT COMBINATION WITH VARIABLE RETARDATION 


We shall now show that a combination of three birefringent plates, 
the first and last of which have their corresponding principal planes parallel 
and their retardations identical, will behave as a single purely birefringent 
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plate of retardation 25, whose fast vibration direction is inclined at an 
angle c, to that of the first plate, where 25 and c, are given by: 


Cos 5 = cos 28, cos 5, — sin 28, sin 8, cos 2c (a) 


and 
cot 2c, = cosec 2c (sin 25, cot 5, + cos 28, cos 2c) (b) 


In the above relations 25, is the common retardation of the first and 
last plates, 25, that of the central plate, and c the angle between the fast 
vibration directions of the central plate and the other two. 


Before proceeding to prove these relations, we may mention that by 
choosing 2 quarter-wave plates as the first and last plates, and a half-wave 
plate as the central plate, we can, by rotating the central plate, vary the 
retardation of the combination continuously from 0 to 27; the principal 
planes of the combination will always be inclined at 45° to those of the first 
plate. 


Referring to Fig. 2, a rotation 25, about AO followed by a rotation 
28, about BO, where A and B give the orientation of the fast axes of the 
first and middle plates respectively, is equivalent to a rotation 2¢ about CO— 
where CAB = 3,, CBA = 8, and 7— ACB = ¢. Let us draw an arc CD 
such that ACD = ¢, D being the point of intersection of this arc with the 
equator. Considering now the triangle ACD, the rotation 2¢ about C 
(representing the combined action of the first 2 plates) followed by a rota- 
tion 25, about A (which now represents the action of the third plate) is equiva- 
lent to a rotation about D through twice the angle external to ADC, Th¢ 
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combination is therefore equivalent to a birefringent plate of retardation 28 
whose fast axis is inclined at c, to that of the first plate—where 5 = 7 — ADC, 
and 2c, = are AD. 


The expressions which we shall derive for 5 and c, do not depend on 
assuming that D lies within AB. 


From the triangle ADC, denoting arc CD by d 

cos 5, = cos ¢cos 5 + sin dsin 5 cos d (1) 
From the triangle DBC 

cos 5, = cos 2¢ cos 6 + sin 2¢ sin 5 cos d (2) 


Eliminating d by multiplying (1) by 2cos ¢ and subtracting from (2), 
we obtain: 


cos 6 = 2 cos ¢cos 5, — cos 5, 
The value of cos¢ to be substituted is obtained from the triangle ABC. 


cos = cos 5, cos 5, — sin 5, sin 5, cos 2c 
Hence 
cos 5 = cos 26, cos 5, — sin 28, sin 5, cos 2c 


which is the required expression (a). 


To get an expression for c,, we use a well-known trigonometrical theorem, 
which gives from triangle ABC, on denoting are AC by b 


cot 2c sin b = cot ACB sin 3, + cos b cos 3, 
and from triangle ACD 


cot 2c, sin b = cot ACD sin 8, + cos bcos 3, 


or 
sin b (cot 2c, — cot 2c) = 2 sin 3, cot ¢ 
Substituting 
sin b = sin 8, and 
we obtain, 
__ 2sin 8, 
sin 2c (cot 2c, — cot 2c) = sin 5, cos 


sin 2c cot 2c, — cos 2c = aa ag (cos 8, cos 5, — sin 8, sin 8, cos 2c) 
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or 
sin 2c cot 2c, = sin 25, cot 5, + cos 25, cos 2c 
which is the required expression (bd). 


4. CONDITIONS FOR ACHROMATISM 


To make a combination of the type described above achromatic in the 
immediate neighbourhood of the wave-length Awe may differentiate the 
expressions on the right-hand side of the equations (a) and (b), and equate to 
zero. Since however our object is to attain achromatism over a finite range of 
wave-lengths we proceed in a different manner. 

Let X’ and A” be two wave-lengths on either side of A for which the retarda- 
tions of each individual plate are multiplied by a factor (1 — f) and (1 + f) 
respectively. Values corresponding to these two wave-lengths will be denoted 
by corresponding single- and double-primed symbols. 


We impose the conditions that the retardation 28’ and 28” of the combi- 
nation for these two wave-lengths be equal to one another and to 24, the 
required retardation of the achromatic combination; also that the orienta- 
tions of the fast axes of the combination for these two wave-lengths should 
be the same, i.e., c,’=c,". We then get three equations involving the 
three unknowns 4,, 5, and c: 


cos 4 = cos 28,’ cos 5,’ — sin 28’, sin 5,’ cos 2c (1) 
cos 4 = cos 28,” cos 5,” — sin 28,” sin 5,” cos 2c (2) 

sin 28,’ cot 5,’ + cos 28,’ cos 2c = sin 28,” cot 5,” + cos 28,” 
cos 2c (3) 


From these three relations we may get two relations involving only the un- 
knowns 5, and 585. 


Thus eliminating cos 2c from (1) and (2) 


cos 8,’ cos 28,’ — cos 4 __ cos 8,” cos 28,”— cos 4 


sin 28,'sin8, 28," sin,” (4) 
For the second equation, we eliminate cos 2c as well as cos 4 between the 


three equations. Thus equating the right-hand expression of (1) and (2), 
and then substituting the value of cos 2c from (3) we get: 


cos 28,’ cos 5,’— cos 28,” cos 8,” __ sin 28,’ cot 5,’— sin 28,” cot 8,” 


sin 28,’ sin 8,'— sin 28,” sin 3,” — cos 28;"—cos 28 
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which on simplification yields, 


(1 — cos 28,’ cos 28,”) (cos 5,’ + cos 5,”) 


= sin 23,’ sin 235” sin 5,’ cos 5,” + cos 5,’ sin 3,” 
. 7 sin 5,’ sin 5,” 
Substituting 
5,’ = (1 — f) 8, and 8,” = (1 + f) 
we get 


sin 28,’ sin 8,” sin 28, 
(1 — cos 28,’ cos 28,”) 2 cos 8, cos f 8, = (5) 
We have to find values of 5, and 5, which simultaneously satisfy (4) 
and (5). The task is considerably simplified when we note that equation (5) 


is identically satisfied for 5, = 7/2. The corresponding value of 5, may be 
found from (4), which now becomes: 


sin f 5 - sin 2 (8,' + 8;") = cos 4 (sin 28," — sin 28,') 
Substituting 


= (1 —f) 8; and 8,” = (1+ f) 
this gives 


sin 
sin f . 28, = - 


(I) 


An expression for cos 2c is got by eliminating 4 from (1) and (2). 
sin f 
2 _ cos 28,’ + cos 28,” 


ma sin 28,’ — sin 28,” 
2 


cos 2c = 


cos f 


tan f° 28, 


cos 2c = 


For any required retardation 24 and for any choice of the arbitrary 
parameter f which determines the range of achromatism, the retardation 25, 
of the first and last plates of the achromatic combination may be found 
from I; since this is a transcendental equation it has to be solved by an itera- 
tion procedure. The orientation c of the central plate, which should be a 
half-wave plate, may then be obtained from II. The orientation c,’ of the 
principal planes of the combination may, if required, be calculated from 


sin 2¢ cot 2ey’ = sin (1 — f) 26 tan + cos(1—f) 28, cos2e 
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The retardation and the orientation of the principal planes of the achro- 
matic retardation plate, taken as a whole, will disperse slightly with wave- 
length. The combination will possess exactly the same retardation 24 and 
the same orientation c,’ of its fast vibration direction, for two wave- 
lengths; these wave-lengths, if we neglect the dispersion of birefringence, 
will be (1 + f) A and (1 — f)A; and we may, somewhat arbitrarily, consider 
the range of achromaiism as extending from (1 + ./2f)A to (1 — V2f)A. 


We shall present the numerical solutions only for the case of an achro- 
matic quarter-wave plate—for which we must substitute 24 = 90°. If we 
choose f = -18 the solution of I can be shown to be 28, = 115° 42’; so that 
from II, cos 2c = — -7639 or 2c = 139° 48’. The corresponding value of 
c;' is got from III, which yields cot 2c,’ = -5489 or 2c,’ = 61° 14’. 


On the other hand if f= -1414 we get 
28, = 115° 30’, 2c = 140° 26’, 2c,’ = 61° 30’ 


Thus an achromatic quarter-wave plate is obtained by superposing three 
plates of the same birefringent material.(the dispersion of the birefringence of 
which need not be negligible). The central plate should be a half-wave plate 
for the wave-length A in the middle of the spectral range to be covered. The 
first and last plates should have their principal planes in parallel orientation 
and their retardations of the same magnitude—equal to, say, 25, for the 
wave-length A. With respect to the common fast vibration direction of the 
first and last plates, let the inclinations of the fast vibration direction of the 
central plate, and of the combination taken as a whole, be c and c,’ respec- 
tively. Then for covering the range of wave-lengths from about 1-25 A to 
75 » the set of values, 25, = 115° 42’, c = 69° 54’ may be used—for which 
cy’ = 30° 37’. The achromatism is not high but can be increased by restrict- 
ing the range to be covered. Thus for covering the region from 1-2 A to 
*8 A the values suitable are: 


28, = 115° c= 13’, 
for which 
Cy = 30° 45’ 


We may mention that the range over which these combinations may be 
considered achromatic would have to be altered if the dispersion of the bire- 
fringence is appreciable, but their essentially achromatic nature would not be 
changed. For retardation plates of mica, the dispersion of bitefringence 
is negligible,* 
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EXPERIMENTAL VERIFICATION 


Two plates of mica were cleaved, whose thicknesses were in the ratio of 
180: 115-75, the former being nearly a half-wave plate for the D lines of 
sodium. The second plate was cut into two portions which were then super- 
posed and cemented together at a corner—the principal planes of the two 
portions being made exactly coincident while the cement was still wet. After 
the cement was dry, the half-wave plate—with a drop of copal varnish on 
either side—was inserted between the other two plates at roughly the calcu- 
lated angle. The combination was placed on a mirror and, with the aid of 
. another mirror, white light from a point source was made to pass normally 
through a polaroid and the combination of plates and then back again— 
after which it reached the eye. The orientation of the central plate was 
then gradually altered while the varnish was still wet—and a position was 
found where the image of the point source could be almost completely 
extinguished by rotating the polaroid. It was also verified that the 
combination thus prepared showed two practically black extinction posi- 
tions when rotated between crossed polaroids. 


The author’s thanks are due to Prof. Sir C. V. Raman for his kind 
interest in this work. 


SUMMARY 


An achromatic quarter-wave plate is obtained by superposing three bire- 
fringent plates of the same material; the first and last should have the same 
retardation 28,, their fast vibration directions being parallel to one another 
but inclined at a specific angle c to that of the central plate—of retardation 7. 
The desired range of achromatism determines the optimum values of 25, 
and c (which, in turn, will determine the orientation of the effective principal 
planes of the combination). As an example, using mica retardation plates 
prepared for Hg 5461, the range from 4100A to 6800A is covered with 
25, = 115° 42’ and c = 69° 54’. 


Further, for a particular wave-length, a birefringent compensator of 
variable retardation (0 to 2) is obtained by interposing a half-wave plate 
that can be rotated in its own plane, between two quarter-wave plates that 
have their fast vibration directions parallel. 


The results follow from the Poincaré sphere by geometrically compound- 
ing successive rotations. 
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INTRODUCTION 


THE subject of inverse probability has unfortunately been regarded with 
general suspicion and so not many attempts have been made to apply the 
theory even to cases where unambiguous answers can be given, on a proper 
mathematical formulation of the problems. It is the object of this paper 
to emphasise that in the case of evolutionary stochastic processes, the applica- 
tion of Bayes’ theorem yields results as meaningful and rigorous as could 


be obtained by the application of the concepts of ordinary (as contrasted 
with inverse) probability theory. 


We shall first formulate the basic problem of inverse probability in 
relation to univariate stochastic processes and point out why it was 
erroneously supposed to be a difficult problem. We shall then show how 
the problem can be solved rigorously by reconsidering the Chapman- 


Kolmogoroff equation and consequently the stochastic differential equations 
derived therefrom, in a new light. 


2. FORMULATION OF THE PROBLEM 


Let us consider a stochastic process representing the ‘ evolution’ in a 
Markovian manner with respect to the one-dimensional parameter ¢ of the 
probability distribution of a continuous stochastic variate x(t). We are 
thus interested in the probability frequency function (p.f.f.) 7 (x,| X93 ty, to) dx4 
representing the probability that x(¢,) lies between x, and x, + dx, given 
that x (to)= Xo (4 > to). 

It shall be assumed that the process is homogeneous in ¢ and so z is a 
function only of the interval t,; — tf) and not t,, f) severally. Though it is 
customary to set fp= 0 and write t, — fy as t,, we shall deliberately refrain 
from doing so, for reasons which will be clear presently. On the contrary 
we shall use the notation 7 (x; ¢) to represent the p.f.f. of x (t) when we do 


* Read at the Annual Meeting of the Indian Academy of Sciences at Belgaum, December, 
1954, 
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not specify the “ initial” conditions, i.e., the p.f.f. at some ¢’ 


Thus by definition we immediately write 


(1) 


When we deal with distributions at two different values of ft, say t, and fp, 
(t,;> to), we shall call the distribution at ¢, the emergent spectrum with respect 
to the injected spectrum at fp, and 7 (x,| X93 t;, fo) the fundamental transition 
probability for the finite interval t, — to. It is quite clear that 7 (x,| Xo; t,, to) 
is only a particular case of 7(x,; 1¢,) where the injection spectrum at fp is 


given by 7 (x; f))= 5 (x — xX») where 4 is the Dirac delta function. 7 satisfies 
the Chapman-Kolmogoroff equation, 


w(x; J 


t= J to) X03 bo) (to< ty). 


to +7) x's ty, to + 7) dx’ (2) 


for all values of 7 such that ty < ty + 7 — ¢t, and in particular 


X03 ty to) m(x'| X93 to + 7, to) ty to + 7) dx’. (3) 


By making ¢) + 7 — #,, we obtain the well-known forward differential equa- 
tion of Kolmogoroff for 7 provided we know the limit of (x,| x’; t, f) 
as t, — t ++ 0, i.e., to speak in physical terms we know what ‘ happens’ 
to the stochastic variable is an infinitesimal interval 4. 

From equation (1) we find that the emergent p.f.f. 7(x,; t,) can be 
computed if we know the injected p.f.f. 7(x9; to) and the fundamental 
transition probability 7(x,| x9; But the computation of t)) 
given 7(x,; ¢,) seems at first a formidable task since it involves the solu- 
tion of equation (1) which is a Volterra integral equation of the first kind. 
It is the principal object of the present paper to show that no complex proce- 
dure is necessary to solve for the injected spectrum in view of the peculiar 
properties of the fundamental transition probability. We shall prove by very 
simple arguments that the inversion of (1) is given by 


Tt (Xo; t= ty) 7 (Xo) X13 thy to) > to), (4) 


where 7 (Xo | X13 t, t) is a function obtained by replacing t, — ty by the 
negative quantity ty — t,, in 7(X| X13 t, t) which we remember is only a 
function of — and not t,, fy severally. unlike has no probability 
significance and is only a functional operator. 


+ Throughout this paper denotes a p.f.f. Distinction between different distributions will 
be obvious from the variable used to define 7. 


: ‘i 
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In equation (1) since (x; and 7(x,| f9) are both probability 
densities, they are positive and so 7 (x; t,) exists for any ft, > fo, i.e., ty — to 
is positive and finite and can be chosen as large as we please. If 
m(x| Xo; t, te) tends to a stationary distribution as t, — ty > + 0©, so 
also does 7(x,; 


On the contrary in (4) since 7 is not a p.f.f. 7(x; f») does not remain 
a p.f.f. for all values of tj << t,. It shall be proved that there exists a tp 
such that 7(x»; ft ) is a p.f.f. in the domain tp < t»<¢t, and the process 
cannot be traced back to a point t < ftp. 


3. SOLUTION OF THE PROBLEM 


It is necessary to recall the method of derivation of the forward differ- 
ential equation for 7(x; t) from the Chapman-Kolmogoroff equation. 
In particular, we first consider the simple case when 


X03 ty, to) > R (x,| Xo) 4 + (x1 — — R (x’| dx’) (5) 
as tt—ty>4 >+0. 


In such a case, (1) assumes the form 


a(x; te + 4)= 4 J (Xo; to) R(x| Xo) 
+ t)) {1 — 4 fR (x’| x) dx’} + 0(4?). (6) 


(6) can be reduced to a differential equation in ¢ by making 4 > +0. From 
(6) + 4) can be computed if 7 (x; f)) is known, and thus the emer- 
gent p.f.f. 7(x,; is obtained from 7(x»; by a process of integration, 
ie., by solving the differential equation in 1. 


We notice we can write (6) in the following form: 
to + Af to + 4) 


+m (X5 bo) — to + 4) + APR x) dx’ + 0(4?) 


(7) 


which can be re-written in the equivalent form (8) by using the device of 
writing t) --.4 and ft, for tf, and t, + 4 respectively. 


to — A= to) + Aa(x; to) J R (x’| x) dx’ 


45 (x95 te) dy @) 


to 
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In this form we immediately realise a(x; ty — 4) can be obtained from 
a(x; to) by merely replacing 4 by — A in w(xy; t+ A). Thus 7 (X93 to) 
can be obtained by a process of integration from m(x,; t,), (t, > to). 


Considering processes more general than that defined by (1), Bartlett 
(see Ramakrishnan, 1952) has observed that the forward differential equation 
for w can be written in the generalised symbolic form, 


(7 positive) omitting for the moment the argument x to indicate that the 
right-hand side may involve several values of the argument. If m(fo) is a 
p.f£.f. (t) + 7) is a p.f.f. for all from 0 to + co. Noting that the above 
equation is derived from the following 


+7 + A=a(ty t+ + 7} -4+0(4%), (10) 
we re-write it in the form 


(to + A= a (to +2) + F fa (to +7 + 4+0(49). (11) 


Therefore using the same device as before of writing ty — 7+ 4 and 
ty — r in the place of t, + 7 and ty +7+ 4, 


a(ty— 7+ 7) = — 4F {a (ty — 7)} 


In other words 7(t, — 7) is obtained from 7 (to) by merely replacing 7 by 
in w(t,+ 7). 


It now remains to show that this process of inversion, i.e., determining 
the p.f.f. at ts — 7 from the p.f.f. at t9 cannot be carried on for all values 
of r except in the trivial case when the emergent distribution at fy is a sta- 
tionary distribution. 


To do this we first consider the simpler case defined by (5) and examine 
the essential difference between (6) and (7). Noting that R (x| x») and 
a(x; to) are non-negative in (6), we find that 7(x; tf) + 4) is also non-nega- 
tive. On the contrary 7(x; to — 4) will be negative for such values of x 
for which m(x; t,)=0 if R(x| x’) is positive for some value of x’. In 
such a case 7 (x; f9) cannot be an ‘ emergent’ distribution, i.e., it cannot be 
traced back. Thus there exists a smallest value ty for r such that (x; te—Tp) 


i 
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is zero for some value of x. Thus tp=to — Tp is the earliest point to which 
we can trace back the process. Identical arguments apply to a 7 function 
satisfying Bartlett’s generalised symbolic equation. We call tp, the ‘ primitive 
origin’ of the process and rp the ‘ age ’ of the process at fy. The distribu- 
tion at any point ty) — r, 0 << 7 < Tp is given by equation (12), i.e., if we 
inject +(x; to — 7) satisfying (12), then it will emerge as 7(x; fo) at fy. 
It also follows that 7(x; ty — 7,), the injection spectrum at t, — 7, is the 
emergent spectrum with respect to 7 (x; ft) — 72) the corresponding injection 
spectrum at fy — 72, (tT, > 74). (x; t — 7) satisfies the modified Chapman- 
Kolmogoroff equation 


a(x; ty — J tg — 7’) w(x| to— 7, to — 7’) dx’, (13) 


where 
ty—m<ty—T< 
> 


The distribution at t) — ty can only be an injected distribution and 
cannot be the emergent distribution of a process started ‘ earlier’. Therefore 
we call it the primitive a priori spectrum. 


Illustration in the case when x(t) is discrete-—We shall now consider 
the case when x(t) is discrete and can assume mutually exclusive values 
Xj, and define 7) as the probability (not probability 
density) that x(t) = xj. Ifw(j|i; ty, to) > dast,— ts > 4 > +0 
it is well known that z satisfies the matrix differential equation [cf. integro- 
differential equation when x(t) is continuous] 


where z (j; ) is a column vector with components corresponding to different 

values of j and [R] is the matrix with non-diagonal elements R (j| i) and the 

diagonal elements defined by R (i|i)= — Y R(j\i). The solution of the 
j 


above equation is given by 
ty) = to ta). (15) 


By considering the infinitesimal transformation 7(j; ¢) to 7(j; t+ 4) we 
note that if, 7(j; t) >O forall j soisw(j; t+ 4). Thus 7(j; t,) > 0 for 
all j for all t, >t provided >O0 for all j. 
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By arguments now familiar from the previous section 7(j; 1 — 4) can be 
obtained from 7(j; t) and so 


(7; (ty > to) (16) 


i.e., given the emergent distribution 7 (j; t,) we compute the injected distri- 
bution at fy. This computation can be done for ¢, in the domain t;—Tp < 
At to=t,— Tp, — at least for one value of j. 


To illustrate this we shall take a system capable of assuming one of four 
states, i.e., j= 1, 2, 3, 4. We shall take the following values for R: 


R(i/1) R({2) -12 0 00 
[na 1) RQ|2) RQ} 3) 


1-8 
R(4\1) R(4\2) 4) 2 


We assume that at some /,, the values of 7(j; ¢,) are as given below: 
then the primitive origin ¢, — tp is determined. 


(2; -1, 7(2; ty — tp)= +2074 

for = -1097 
t= -2, #3; — 0 ( 
(4; t)= 7(4; — -4197 


Conditional Inverse Probabilities—We have observed that the p.f.f. 
(Xy| Xo3 ty fo), > to) is a particular case of m(x,; with a(x; to) 
= 8(x — x,). In this case the process cannot be traced back further than 
to, i.€., to a point t < fy, since to) = 0 for x# Xo. Hence w(x; ty — 4) 
is negative (4 positive) for xAx,. Therefore a delta function can only 
be an injected spectrum and cannot be the emergent spectrum of a stochastic 
process. Thus 7 (x| to), (¢ > has no probability interpretation. 

On the other hand, in a previous paper, the author (1954) defined P 
(X9| X13 as the conditional inverse probability that x lies between 
and xX9+dx, given that x (t,)=,, (to < t,). This has to be distinguished from 
z. When in defining P we state that x (t,) = X,, we mean that the observed 
value of the emergent spectrum at t, is x, and not that the emergent spectrum is 
a delta function. In we require the emergent spectrum to be a delta function 
and this is not possible. Hence z is not a probability magnitude. 
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_In his previous contribution, the author showed that P (x9| x1; to, 4) 
cannot be obtained unless the injected spectrum at a certain ¢ <fy is known. 
If that is known to be m(x; f) then by considering the joint probability of 
obtaining x (t)= x, x (to) = Xo and x (t,)= x, and applying Bayes’ theorem, 
it was shown that 


f(x; t) x3 4, 1) dx 


(X03 to) X03 Fo) 

= &) t<t%<t) (1%) 
Ifr(x; t) is known 7 (Xp; fo) and 7(x,; ¢,) can be computed from the funda- 
mental transition probabilities corresponding to the intervals t, — ¢ and 


t;— tp. Hence the method of computation of P described in that paper 
is equivalent to the present one. 


In this paper we have shown that if 7(x; ft) is known at some point, 
it can be extended to the entire permissible domain of t and so P can be 
computed if the distribution at some ¢, and of course the fundamental transi- 
tidn probability for a finite interval are known. 


PROBLEMS INVOLVING A PROBABILITY FREQUENCY FUNCTION FOR 
THE PRIMITIVE ORIGIN 


Till now we have assumed that the distribution at ¢ is emergent with 
respect to that at t — 4 (except when ¢ is the primitive origin) and injected 
with respect to that at ¢-++ 4. Or, in other words, given the distribution at 
to, the injected distribution at ¢ < fy is a function of ty, — ¢. The primitive 
origin therefore is a definite point and is uniquely determined if we know 
the probability frequency function at fy. 


We shall now consider the following class of problems. We are 
given the probability frequency function at ft). The primitive origin 
of the process lies between /,— 7 and ty —(7-}+ dr) with probability 
n(7)dr. At this origin a p.f.f. ¢(x) independent of 7 is injected. Given 
»(r) our problem is to determine the injection spectrum ¢ (x). 


Using simple probability arguments we write 
Pdr (x) x5 toy 7) (2) (18) 
0 


Since f)) and x’, to, 7) are known, or @(x’) can be deter- 
mined if the other is known, by the classical theory of integral equations; 
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but the explicit solution of ¢(x’) or »(t) amenable to numerical computa- 
tion is a very difficult problem. On the contrary, the computation of 7 
given ¢ is a much easier problem. Considering the very simple special case 
when x is a discrete variable representing a Poisson distribution 
(x| x 3 & t— Tt) x’! 
and is independent of 1, 


w(x; t= fer a(x) det (19 


As expected 7(X; fo) is independent of f, since (x| x’; t, t — 7) is only a 
‘function of 7. The computation of 7 given 7 is still a difficult problem. 
Consael has discussed a class of distributions defined by 


The physical meaning of this class of distributions is obvious. The 
parameter A of the Poisson distribution corresponding to an interval t is 
itself a stochastic variable with a p.f.f. y(A). Comparing (20) with (19), 
we find that all the results of Consael apply to processes defined by (19) 
provided we replace A by ¢, since the Poisson distribution is symmetrical 
with respect to ¢ and A. 


Problems involving a “‘switch-off” of the process.—Let ty be the primitive 
origin of a process where we inject the primitive p.f.f. 6(x). Let us assume 
that we ‘‘switch-off” the process between tj, + 7 and t, + 7+ dr with 
probability » (7) dr. After “* switch-off”’ the stochastic variable retains the 
value it had at r till infinity. Thus the p.f-f. at tg + 7 as 7 + 00 is given by 


m dx $ (x) (x'| x3 to + 7, te) dt (21) 


Comparing this equation with (18) we find z (x’ ; 00) is identical with = (x; 
to) of the previous problem. It is interesting to note that in the previous case 
the point of origin has a p.f.f. while here the point of switch-off has a p.ff. 
and therefore the duration of the process has the same p.f.f. in both cases ! 

Incidentally the author wishes to mention that the process of acceleration 
of cosmic ray particles in the galactic magnetic field proposed by Fermi is 


a process which can be interpreted to belong to either of the two equivalent 
types discussed above. 


¢ We assume, without loss of generality, that x’ = 0. 


re 
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Process involving absorption”’’.—Finally, we advert to an interesting 
problem in inverse probability occurring in physics (relating to the range 
of a fast particle passing through matter) which for some unknown reason 
has not been investigated in a rigorous manner. 

A fast particle of initial energy E, passes through matter and loses 
energy through ionisation or radiation. When it reaches an energy zero or 
some critical value E, we say the particle is ‘ absorbed’ or stopped, and the 
distance ¢ of penetration (considering the passage to be one-dimensional) 
till it is absorbed is called its range. The author has, in a separate contribu- 
tion to this journal, derived an expression for P(E,; ¢) dt, the probability 
that a particle of initial energy E, is absorbed between ¢ and ¢ + dt if we 
know the fundamental transition probability 7(E|E,; ¢, 0). 

If we now assume that the initial energy of the particle is a stochastic 
variable defined by the p.f.f. ¢(Eo), i.e., if (Eo) is the primitive injected 
spectrum, the p.f.f. of the range is given by 


P()dt= P(E, ; 2) dB, (22) 


Given P (t) and P (E,; 1) it is theoretically possible to compute the primi- 
tive injected spectrum ¢(E,) but this is not as easy as the computation of 
¢(E,) if we are given the emergent spectrum at f. 

There still remains the further question of conditional inverse probability. 
Given that the observed range was ft, to compute the probability ¢4;, (Ey) dE, 
that the initial energy of the particle lies between E, and E, -++ dE. By the 
simple application of Bayes’ theorem we find ¢;, (Eo) is given by 


$ (Eo) P(E; t)__ P (Eo; #) (23) 


Ed = PE; 0) dE, 

It is clear that the conditional inverse probability cannot be obtained 
unless we know P(t) and P(E); ¢) or @(E,) and P(E,; f). In the opinion 
of the author, the function ¢;, (Eo) is physically more important than P (Ep; f) 
since in many experiments the energy of the particle is estimated from the. 
range and somewhow this aspect seems to have attracted little attention till 


now. 
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ABSTRACT 


Small local dislocations in processed emulsions give rise to spurious 
scattering whose effect on tracks is indistinguishable from the multiple 
scattering produced by coulomb interaction. Because of this effect 
the atomic number of fast primary cosmic ray nuclei will be underesti- 
mated in more than 50% of the cases in experiments where multiple 
scattering of tracks and ionization measurements are employed for charge 
determination. The spurious scattering was investigated by measuring 
the track contours of 100 very long tracks due to energetic primary 
particles; the effect was present in comparable strength in all plates and 
all types of emulsions which have been investigated. It is too small to 
affect measurements on medium energy particle tracks (protons with 
energy below 600 MeV and heavy nuclei with energy below 300 MeV/ 
nucleon). It dominates, however, other sources or error and noise for 
tracks of energetic particles although it does not preclude occasional 
observations of very low scattering values. New methods for measuring 
various forms of noise have also been developed in the course of this 


work and the noise level for scattering measurements has been reduced 
below previously accepted values. 


Spurious scattering is presumably largely responsible for the discre- 
pancies which appear when one compares the primary charge and energy 
spectra derived from experiments involving scattering measurements 


with the corresponding spectra derived from numerous other experiments 
which employ different techniques. 


If the experiments based on scattering measurements are omitted, 
the remaining evidence strongly favours a spectrum in which the energy 
per nucleon is nearly independent of atomic number for all primaries. 
It also favours a charge distribution which has a pronounced minimum 
for charges 3< Z< 5 and, therefore, yields a fairly low upper limit for 
the amount of interstellar matter traversed by primary cosmic ray nuclei. 


1. INTRODUCTION 


THE determination of multiple coulomb scattering of particle tracks in 
nuclear emulsion has been widely used in conjunction with ionization or 


| 
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range measurements for obtaining information on mass, charge and energy 
of particles in cosmic ray experiments. The conclusions derived from such 
scattering measurements can often be confirmed by independent evidence 
from other sources. At other times, however, conclusions derived from 
scattering measurements differ from those obtained by other techniques. 
Such a discrepancy arose for instance, in experiments designed to determine 
the charge and energy distribution of heavy nuclei in the primary cosmic 
radiation. 

Experiments not involving scattering measurements have been carried out 
with nuclear emulsions,’ ionization chambers,® proportional counters,® 
cloud chambers triggered by proportional counters,’ and Czerenkov detectors 
incorporated into an ordinary counter telescope*; they support or are con- 
sistent with the conclusion that all components of the primary cosmic radia- 
tion have closely similar energy spectra (when expressed ‘in terms of energy 
per nucleon) and that next to protons and a-particles the group of elements 
comprising Carbon, Nitrogen and Oxygen is strongest while Lithium, Beryl- 
lium and Boron nuclei are either rare or absent. 


On the other hand experiments involving measurements of multiple 
scattering in nuclear emulsions lead to an energy spectrum in which the 
heavy primary component falls off more rapidly with energy than the total 
primary radiation and yield a charge spectrum in which Li, Be, B nuclei are 
approximately as abundant as C, N, O nuclei.*-!? 


This discrepancy has led us to re-examine the limitations of the existing 
techniques for the measurements of multiple coulomb scattering in nuclear 
emulsions. We found that the discrepancy is apparently due to an imperfec- 
tion in the nuclear emulsions now in use. This result seems at first sight 
surprising in view of the fact that numerous experiments have established 
the reliability of conclusions drawn from scattering measurements on tracks 
in nuclear emulsions. In particular it has been shown that the mass of slow 
particles can be correctly derived from the multiple coulomb scatiering of 
their tracks in conjunction with measurements of range or specific ionization. 
It has also been shown that scattering measurements give the correct answer 
when applied to tracks of fairly fast particles of known mass and energy, 
accelerated in the laboratory. Finally, it has been shown that occasionally 
tracks of particles which are believed to initiate high energy nuclear interac- 
tions exhibit very small multiple coulomb scattering, so small that one can 
set the lower limit for the energy of these particles at 10 BeV or even higher, 


However, the imperfection which we have found to be the property of 
all emulsions so far investigated and which gives rise to Spurious Scattering, 


if 
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does not invalidate any of these well established conclusions. The spurious 
scattering is too small to interfere with measurements on slow particles and 
sometimes its effect becomes small for very flat tracks. It does, however, 
play an important role when measuring tracks of particles whose energy 
exceeds 0:6 BeV (pf = | Bev/c) in case of protons and 0-3 BeV/nucleon 
(p8 = 0-5 BeV/c) in case of heavy ions, especially if these particles enter the 
emulsion with dip angles larger than 2°. 


We have found this spurious scattering effect present in all emulsions 
so far investigated: in glass backed emulsions, varying in thickness from 
250 to 1,000, and in stripped emulsions belonging to various batches. 
We have also found it in plates developed in various laboratories (Bombay, 
Bristol, Chicago and Melbourne) and in emulsions with measurable distor- 
tion, as well as in others in which distortion was extremely small. There 
seems to be no correlation between general distortion and spurious scatter- 
ing. So far we have restricted our investigations to Ilford G-5 emulsions. 
We have not yet come across plates which are free from the spurious scatter- 
ing effect or in which the effect is substantially smaller than average. 


The small-scale dislocations which give rise to spurious scattering are 
illustrated in Fig. 1. This figure shows the contours of two a-particle tracks 
which run closely parallel and are never separated by more than 40 py. (The 
ordinate scale in Fig. | has been chosen so as to amplify the y-co-ordinate 
4,500 times.) The two a-particles belong to a shower containing six rela- 
tivistic a-particles which survived the destruction by collision of an energetic 
heavy primary nucleus. The contour of each track was measured in 250» 
steps, but the tracks were measured on different parts of the microscope 
stage motion. The very close correspondence between these two com- 
plicated contours can, therefore, not be attributed to irregularities in the 
screw which moves the stage. It can obviously not be attributed to causes 
which affect each track independently like reading errors, grain noise or 
coulomb scattering of the particles. It is an illustration of the extremely 
complicated small-scale dislocations which occur in emulsions, quite different 
from the C or S-shaped distortions which one meets occasionally. The 
effect illustrated in Fig. 1 exists in all emulsions and for all tracks which we 
have investigated so far, although in the example which we have chosen for 
illustration it is somewhat more pronounced than average. Still most of the 
dislocations common to both the tracks are not larger than fractions of 
microns, and the overall S or C-shaped distortion is moderate (50 covans). 


It is the purpose of this paper to study the effect of these dislocations 
and of the resulting spurious scattering on experiments designed to measure 
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multiple coulomb scattering of particles. This is done in Section III. As 
far as we know, this effect has not previously been discussed and has not been 
considered in the numerous papers which deal with the multiple scattering 
technique in nuclear emulsions. 
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Fic. 1. The contours of two parallel tracks separated by ~ 40” were measured at 250u 
intervals on different parts of the microscope stage. The tracks belong to two (out of six) 
relativistic a-particles and are fragments belonging to a heavy nucleus with 15-20 BeV/nucleon 
destroyed in a nuclear collision. Scattering data are givenin Table IV. Correlation between 
the two contours is due to the local dislocations which give rise to spurious scattering. The 
ordinate in Fig. 1 has been amplified 4,500 times. 


In order to study the effect of spurious scattering we have discussed 
in Section II, various other limitations on measuring accuracy which arise 
in the form of ‘‘ noise” and are fairly well known to workers in this field. 
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We believe, however, that we have made some improvements in the methods 
of measuring these various forms of noise and achieved an appreciable 
reduction in noise level below that previously reported in the literature. 


In the last section (Section IV) we consider the effect of spurious scattering 
on the group of experiments designed to measure the energy and charge 
distribution of heavy primary cosmic rays. We briefly review the existing 
evidence and recalculate the upper limit of interstellar material traversed 


by the cosmic ray beam in the light of present knowledge of its chemical 
composition. 


II. Noise CORRECTION IN SCATTERING MEASUREMENTS 


Throughout this work we have used the co-ordinate method of scattering 
developed by Fowler.'* A track is divided into equal cells and measurements 
consist in determining the distances y, of the dividing points from a straight 
line which is reasonably parallel to the track. The mean absolute value 
of the second co-ordinate difference D is defined by: 


D = </ yn — + Ynie| > (1) 


where < > designates an average over all readings. D contains informa- 
tion on multiple coulomb scattering but must be corrected for various 
errors arising from optical and mechanical imperfections of the scattering 
medium and measuring equipment as well as from errors of the observer. 


We shall discuss the errors in the following order :— 


(a) The temperature noise «; which is caused by differential thermal 
expansion of microscope parts during the measuring period. 


(b) The reading noise €g which arises from the failure of the observer 
to place the goniometer hair-line accurately through the centre 
of gravity of the row of grains or the solid column of developed 
silver which defines the track in the field of view. 


(c) The stage noise «, which arises from imperfection in the motion of 
the microscope stage, and focussing control. 


(d) The grain noise «g. This can arise from the fact that, because the 
grains have a finite size, the centre of gravity of a finite number 
of grains in a row will not in general coincide with the particle 
trajectory. In the case of tracks consisting of solid columns of 
developed silver there may be other effects which cause local devia- 
tions between the centre of the silver column and the particle 
trajectory. One effect may be contagious development which 
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could lead to irregular patches, where the track is broadened. 
Another source of such noise is perhaps due to irregular shrinking 
of the emulsion in the presence of the large inhomogeneity presented 
by the track itself. We shall designate by grain noise any deviation 
between the centre of gravity of the track and the trajectory which 
varies in distances small compared to about 20p. 


(a) Temperature Noise-—The effect of temperature noise can best be 
illustrated by considering Fig. 2. The graph shows the co-ordinates 
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DIP ANGLE core O = e7 
1 ATOMIC NUMBER Z = 12 
1 Cp MEAN READING ERROR 0-025 m4 
1 €, MEAN READING NOISE 0063 
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DISTORTION 30 COVANS 
' 
t 
3 | 
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CELL SCATTERING 
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Fic, 2. The track contour due to a magnesium nucleus * 20 BeV/nucleon} measured 
twice in succession proceeding in the direction indicated by the arrows. The ordinate has - 
been amplified 9,000 times. 


measured at 250 intervals of a track in a 600 Ilford G-5 emulsion. The 
nature and energy of the particle which produced the track and other details 
are not relevant at this point; they are given in the figure. In order to 
exhibit the detailed structure of this fairly straight track the ordinate has 
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been enlarged 9,000 times. The lower curve has been obtained by moving 
along the track in the direction indicated by the arrow; the upper curve 
represents the results of repeating the measurements of the co-ordinates at 


the same points in reverse order, immediately after the first series of measure- 
ments had been completed. 


The two curves shown in Fig. 2 do not coincide but show a divergence 
growing with the time interval which elapsed between the corresponding 
readings. This shift is due to a differential expansion between the arm and 
the stage of the Koristka microscope. The case illustrated in Fig. 2 cor- 
responds to a differential temperature change of about 0-2°C./hour. This 
condition of temperature stability can be attained in an airconditioned room 
where all lamps are far removed from the microscope and the microscope 
illumination (with water filter) has been switched on several hours before 
taking readings. Figure 3 shows the relative motion of optical axis and 
stage when a 60-watt table light is switched on and later switched off at a 
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Fic. 3. Motions of the optical axis of the microscope due to changing temperatyr¢ 
gradients, 
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distance of 12” from the microscope. The differential temperature change 
right after switching on or off is 12° C./hour and even 15 minutes later it is 
still greater than -75° C./hour. 


In a Koristka microscope the temperature noise is approximately given 


by: 


éy = 3-8 At microns (2) 


where T is the temperature difference between microscope stage and arm in 
centigrade, ¢ is the time, and Ar is the mean variation in time intervals between 
successive readings. Thus, even in a room which is not airconditioned, 
and where temperature variations between different parts of the microscope 
of the order of 2° C./hour occur, it is possible to make the temperature noise 
ey negligible by taking successive readings at regular intervals, but the regu- 
larity requirements are then fairly stringent. (Thus, for tracks of optimum 
thickness where as shown below, one can reduce noise from other sources 
to 0-058 », the contribution from thermal effects will dominate unless suc- 
cessive reading intervals differ on the average, by less than 25 seconds.) 


In any case scattering measurements should not be started before the 
microscope has reached thermal equilibrium with its surroundings. In this 
investigation we made our measurements with temperature drifts of the order 
of 0-15° C. per hour and took readings at regular intervals differing by less 
than 15 seconds; under these conditions the effect of temperature drift on 
second differences can then be neglected. 


(6) Determination of Reading Noise ¢,.—Returning again to Fig. 2, 
one sees that while the two track contours are very similar, the co-ordinates 
of a particular point on the track, measured at two different times, differ 
slightly, even when we allow for the temperature drift which has occurred 
during the time interval. This difference between measurements taken at 
the same point along the track and at the same setting of the stage, is 
obviously independent of the shape of track and of the stage motion; it can 
only arise from errors in reading the track co-ordinates. It is easy to isolate 
and measure this source of error even in the presence of large temperature 
drifts. 

Let ay and by designate co-ordinate values obtained for the same points 
on the track in the two different sets of measurements. We can subtract 
corresponding readings and form the second differences in order to avoid 
temperature effects. We then obtain for the reading noise ¢, (which is the 
mean absolute second difference of individual reading errors «,): 
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N 
| (an — + — (bn — + base) | 
= V2 € 


(3) 
Thus, ¢€g is the contribution of individual reading errors to the mean absolute 
value of the second difference D. 


More generally one can obtain from two sets of measurements of the 
same track the reading noise not only for second differences, but in all orders 


In terms of the mean reading error ex of individual co-ordinate measurements. 
this noise is given by: 


& = < | an — ba | > = V2 We, 


(4a) 
8, = < | (Qn — bn) — — | > = V4 MR (4 
= < | (an — bn) — 2 — Snir) + — | > 
= = V2 (46) 
83 = < |(an — by) — 3 — bnis) + 3 (nse — Onis) 
— (dnis — bnis) | > = V40 ex (4d) 
etc. * 


Experimentally we should expect the zero order noise to be strongly 
influenced by temperature changes which occurred in the time interval 
between the two sets of measurements. First order differences should only 
be affected if the temperature changes are large and second order differences 
only if the temperature changes are non-linear. It is interesting to compare 
the theoretical noise ratios in various orders with the experimental ratios. 
The close agreement between the theoretical and experimental values is shown 
in Table I and proves not only that as expected the reading errors “ €, have 


* The numerical factors are easily derived from the integral below: Let x; be variables 


with a probability distribution P (x,;) = exp and let a; be constants. 
V <7 Xp 
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a gaussian distribution but also that in our experiment temperature drifts 
have only a small effect on first order differences and no effect on the deter- 
mination of reading noise in second or third order. 


TABLE I 
| | | 
1¢ 35 
= +577 =1-825 =3-415 
Experimental ..| °659-+:-04 1-85-+-11 3-51-20 


For each track we have determined the reading noise by measuring the 
contour twice and forming second and higher order differences of the discre- 
pancies between the two series of measurements. The experimental values 
obtained in this way depend of course on the thickness of the track which is 
being measured. The dashed curve in Fig. 4 shows the experimentally 
obtained values of ex, for relativistic nuclei of various atomic numbers in 
Ilford G-5 emulsions. The lowest value «xg = 0:042 py is reached for atomic 
numbers Z = 4 or 5; for lower and higher atomic number the reading noise 
increases. It is interesting to observe that the minimum noise corresponds 
to a mean reading error ‘eg (the accuracy with which the centre of the track 
can be located) of only 170 Angstroms. 


(c) The Stage Noise <«s.—When'the stage of the microscope is moved 
with or without simultaneous motion of the focusing fine control, the optical 
axis of the microscope traces a line on the emulsion surface. If this line is 
not straight, the motion of the microscope controls will give rise to stage 
noise. In order to measure this stage noise we have adopted the following 
procedure. The contour of a track is measured over a distance of 9 mm. 
using successively different sections (i) of the screw which produces the stage 
motion. The contour over each section is measured twice in order to deter- 
mine the reading noise «g by the method described above. We chose a 
track due to a relativistic particle of charge Z = 4 for which the reading 
noise €z was found to be exceptionally low. We then have the following 
relations :— 


(5) 
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Fic. 4. The reading noise ¢, and the combined reading and grain noise €,,, in Ilford 
G-5 emulsions for tracks of relativistic particles with different atomic numbers. 


where the first relation represents second differences obtained from two 
series of measurements, carried out on the same section of the screw, but 
at different times, while the second relation refers to second differences 
obtained from measurements carried out over two different sections of the 


screw. €x,s5 represents the combined noise due to reading error and stage 
motion. 


We find that in our Koristka microscope ex is equal to €,,; within the 
statistical error for all values of cell lengths upto 2,000, and has for the 
particular track chosen for these measurements the numerical value 0-037 p. 
We conclude that the stage noise increases our reading errors by less than 
0-004» or 40A, and is negligible in all scattering measurements. 


In Fig. 5, the combined reading and stage noise ¢,,, and the total noise 
€x,c,s (including grain noise) for fast particles of charge Z = 1 and Z=4 
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Fic. 5. The combined reading and stage noise and the total noise for tracks due to fast 
particles of charge Z=1 and Z=4 is measured on our Koristka microscope for various cell 
lengths and is compared with corresponding values of singly charged particle tracks obtained 
by Gottstein! on microscopes in other laboratories. 


measured on our Koristka microscope is compared with corresponding 
values for singly charged particles obtained by Gottstein“ on various other 
microscopes in use at’ Bristol, Géttingen and Milano. 


(d) The Grain Noise ¢«g.—If the tracks whose scattering is to be 
measured have many gaps, the position of the trajectory has to be estimated 
from the centre of gravity of the small number of grains which define the, 
trajectory in the field of view. Since the centres of individual grains will be 
distributed in a random fashion close to the trajectory, we must expect the 
resulting noise to be inversely proportional to the square-root of the grain 
density (eg = A/+/n). Such grain noise is, therefore, much more important 
for scattering measurements on fast protons than for measurements on 
heavy nuclei. However, there are other effects which can produce localized 
deviations between the particle trajectory and the centre of gravity of the 
track, like irregular broadening by contagious development or irregular 
shrinkage of the emulsion in the presence of the inhomogeneity produced 
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by the track itself. Such effects may increase with ionization. We define, 
therefore, grain noise ¢«g as the noise due to any irregularity in the distribu- 
tion of developed silver around the particle trajectory provided that these 
irregularities occur in distances short compared to 20». 


The following method was used for determining the grain noise. We 
compare again the contours of the same track obtained in two different series 
of measurements; however, this time the points at which the contour is 
determined in the second measurement series are displaced along the track 
by 20, with respect to the points used in the first series. Thereby, a differ- 
ent set of grains is used to determine the position of the trajectory in the 
two series of measurements. If we designate displacement by an asterisk 


we have 
D — D* = V2 €R,G (7) 


In this way the grain noise is slightly overestimated because spurious 
or real scattering along the 20, displacements may also contribute to the 
difference between the two contours. But this error is negligible except 
for tracks whose spurious or real scattering is considerably higher than the 
tracks investigated here (pB < 100 MeV/c.). 


The reading noise €y, the combined reading and grain noise €g ¢ and 
the grain noise (calculated from the relation =€g?+ €,*) are plotted 
in Fig. 4 for relativistic particles of various charges. The steep decline of 
«, between Z = 1 and Z = 22 is certainly an effect of the increased grain 
density and the resulting decrease in the mean distance between the centre 
of gravity of the track and the trajectory. The gradual rise of «, for higher 
values of Z must be attributed to one of the other causes mentioned above. 


One may summarize the discussion in this section as follows :— 


(1) The effect of relative thermal motions of various parts of the micro- 
scope during the period in which measurements are carried out is not in 


general negligible but has been rendered negligible in our experiment by the 
following precautions :— 


(a) The measurements were carried out in a temperature controlled room; 


(b) All heat sources were removed from the neighbourhood of the 
microscope ; 


(c) The microscope was permitted to reach thermal equilibrium after 
the illumination had been switched on; and 


(d) Successive readings were taken at regular intervals. 


(2) Stage noise in our Koristka microscope is too small to affect 
measurements even for those tracks which permit the optimum accuracy. 
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Under these conditions we can write: 


D= + 4?)t (8) 


where as before: 
D is the mean absolute value of the second co-ordinate differences ; 
€x,¢ is the combined reading and grain noise; and 


A contains the contribution of spurious and real scattering to the 
measured second differences. 


For singly and doubly charged particles, and for particles with atomic 
number Z> 8 where the grain noise is large, we always repeat the contour 
measurements using a displacement of 20» as explained before. One obtains 
then from Equations 7 and 8: 


A? = (D)? — 4 D—D*)? mac) 


For tracks due to particles of charge 2< Z< 8 we sometimes used 
the same procedure but more often we repeated the contour measurements 
without displacement and used for ¢g the constant value 0-038, which as 
one can see from Fig. 4, does not introduce measurable errors. This latter 
procedure leads to: 


A? = (D)? — $(D — D’)? — (0-038)? (10) 


Equations 7, 9 and 10 are equally valid if one uses third or higher order 
differences rather than second differences, provided the constant 0-038 is 
multiplied by the appropriate factors (/20/6, 1/70/6, etc.). The use of 
third differences eliminates from scattering measurements most of the effects 
of ordinary large-scale emulsion distortions.” 


EXPERIMENTAL PROCEDURE AND RESULTS 


The following types of tracks were selected for measurement: 


Multiply (and a few singly) charged particles whose energy from 
subsequent nuclear interactions was known to exceed 5 BeV/nucleon; 


Heavy primary nuclei incident at geomagnetic latitudes A< 19°. 


Since it has been shown previously* that all heavy primary particles 
arriving at a given latitude have energies in excess of the geomagnetic cut- 
off energy (the energy required for bare nuclei approaching the earth’s orbit) 
both criteria exclude particles below 5 BeV/nucleon. 
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Sixty tracks have been obtained from a block of 600 » stripped emulsions 
flown from Delhi in 1954. The remaining 40 tracks came from many differ- 
ent types of nuclear emulsions. These include two additional blocks of 
stripped emulsions flown and processed at different times, and glass backed 
emulsions with emulsion thicknesses varying from 250, to 1,000 p. Among 
the latter were plates kindly supplied to us some time ago by Dr. Camerini, 
and others by Dr. Anand, which had been processed at the H.H. Wills 
Physical Laboratory at Bristol, as well as plates processed at the University 
of Chicago and at the University of Melbourne, for which we want to 
express our thanks to Professor M. Schein and Dr. V. D. Hopper, respec- 
tively. 


The geometric selection criteria were based, not on the length of the 
tracks in the emulsion, but on their angle of dip. Comparable numbers of 
tracks were selected with dip angle values in the four intervals: 


8 < cotg@< 12 
12 < cotg @< 20 
< 


20 < cotg @ < 30 


30 < cotg @ 


The contour of each track was measured twice in succession using 250» 
steps. If the tracks were due to singly or doubly charged particles or to a 
particle of charge Z > 8, the two series of measurements were shifted by 
20, with respect to each other and the scattering was determined from 
Equation 9. For tracks of other particles both series of measurements were 
in most cases carried out with the same setting and the scattering deter- 
mined from Equation 10. The scattering value was calculated using second 
and third order co-ordinate differences. In some cases third order differences 
gave lower scattering than second order differences, probably because of 
S or C-shaped distortion. In those cases the value was calculated from 
third differences. These cases were mostly confined to the groups with 
dip angle cotg@ < 20. 


Since all tracks investigated here are very long, statistical fluctuations 
are of no great importance in the results presented below. The consistency 
of measurements on individual tracks can best be judged from results 
obtained on the same track by repetition of measurement. An example 
chosen at random is given in Table II. 


= 
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TABLE II 


Results of Independent Scattering Measurements on a Track due to 
a Particle of Charge Z = 12 


\ Cell D-values (1) 
\ length 
250 500 750 1,000 
1 0-0805 0- 1300 0-1875 0-2325 
2 0-0845 0-1265 0-1910 0-2320 


In Fig. 6 we have averaged separately for each group of dip angles the 
mean scattering value 4 and plotted it against cell length. The same 
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Fic. 6, The second orcer co-ordinate differences due to spurious scattering are plotted 
as a function of cell lengths. Results have been averaged for all tracks belonging to the indicated 
interval of dip angles. The dotted curves show the co-ordinate differences expected from true 
multiple coulomb scattering of heavy nuclei with energies 0-7 and 5 BeV/nucleon respectively, 
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graph shows the expected true scattering value 4; for heavy nuclei of energy 
0-7 and 5 BeV/nucleon. 


From Fig. 4 and Fig. 6 one can see that the average spurious scattering 
greatly exceeds noise values and that spurious scattering increases with cell 
length as rapidly as true scattering. It becomes also apparent that on the 
average flatter tracks have somewhat lower spurious scattering than steeper 
tracks in the interval of dip angle | -9-7°, but below 1-9° we find that spurious 
scattering becomes fairly independent of dip angle. 


Although Fig. 6 exhibits the general trend of spurious scattering with 
track inclination and with cell length, it does not yet convey an idea of the 
spread of values on individual tracks. This is shown in Fig. 7a and b. 
The abscissa in these figures is the experimental scattering value a’) as 
marked on top of the graph and the corresponding energy per nucleon for 
heavy primaries marked on the bottom.7 The ordinate shows the distribu- 


tion in dip angles. (All tracks with cotg@> 65 have been plotted at 
cotg @ = 65.) 


The graphs are divided into three regions marked Z, Z—1 and 
Z-—2 respectively. These regions have the following significance: 


If we had no a priori knowledge that the track under investigation 
was caused by a relativistic particle of energy above 5 BeV/nucleon; 
and 


if we, therefore, attributed the measured scattering to true multiple 
coulomb scattering; and 


if we tried to determine the atomic number of the corresponding 
nucleus from 4-ray density and scattering; 


Then we would obtain the correct charge value for those tracks which 
fall into the region Z; the charge of those which fall into regions 
Z — 1 and Z — 2 would be underestimated by one and two units, 
respectively. 


Figure 7 a shows the distribution of scattering results using 2504 and 
500 » cells; Fig. 7b shows the results obtained with 750 and 1,000» cells. 


On the abscissa we have indicated the limits set to the measurements 
by the various forms of noise discussed in Section II, as well as the energy 
limit imposed by the geomagnetic cut-off. In the absence of spurious 


+ We have used the scattering constant K given by Gottstein etal. 1® Since we are 
dealing with spurious, not real, scattering we have not used any farge angle cut-off, but the 
effect of large angle deflection is in any case unimportant. 
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Fic. 7. Each point gives the dip angle, the spurious scattering parameter a’ 19) and the “energy” 
derived from it, for one of the tracks investigated. Different symbols refer to different cell 
lengths employed. The abscissa also shows the highest energy measurable with 250% and 
500 » cells in the presence of noise. The regions marked Z, Z-1, Z-2 contain the tracks for 
which the atomic numbers derived from ionization plus scattering measurements is underestimated 
by zero, one or two units respectively. The true multiple coulomb scattering for all 100 tracks 
should place each point to the right of the line marked geomagnetic cut-off. 


scattering the points in Fig. 7 a should all lie on the line determined by the 
energy limit set by noise and in Fig. 7 b, all points should lie to the right of 
the energy limit imposed by the geomagnetic cut-off. Almost all points, 
however, lie far to the left of these limits. 
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Together Fig. 7a and b show that spurious scattering will dominate 
completely over other sources of noise for tracks due to particles of a few 
hundred MeV/nucleor. The effect is somewhat more pronounced for 
steeper tracks and shorter cell lengths but, for dip angles below 1-9 
(cotg 6 = 30) and cell lengths above 500,, it becomes nearly independent 
of these variables. 


The graph furthermore indicates that scattering values corresponding 
to energies of 7 BeV/nucleon for heavy nuclei or 14 BeV for protons are 
occasion illy obtainable with cell lengths as low as 1,000. 


In Table III we have listed details of measurements on those tracks 
which produce nuclear collisions indicating energies in excess of 50 BeV/ 
nucleon, that is ten times higher than the minimum value accepted in this 


TABLE III 


Spurious Scattering of Tracks of Particles with Energies 
Above 50 BeV/Nucleon 


Spurious Scattering 
. True Energy | Reading+ for cell lengths 
Description of the Atomic | - Grain Noi ing J; 
late No. Z gy 
BeV /nucleon| éxg (u) 
2504 | 5004 | 10004 10004 
| | 
| | | | 
N-stack 2000 -058 ‘121 +264 -00087 


(600 stripped) 


Melbourne l 1000 -088 | <-088 169-099 00149 
(400 u g.b.) | | 

-058 | -118| -416| 1-370}  -0025 

N-stack 2 300 072 | <-040| -100) +220, -0025 

(600 stripped) | -058| -073| +135 | -0025 


Melbourne 150 -081 | -O71| +170 +388) -0100 
(400 uw g.b.) | | 

077 | «+107 | +324) +775 -0088 

059 | | «+213. +0088 

| <-040| -137| 0083 

500 stripped) -068 <-040 | < -040 0083 

-061 -062 -133| +329 +0083 

‘051 <-040 | 070 | -408 -0088 

N-stack | | -055| -141] -304| -o197 

(600 stripped) +139 | +360} -0127 

N-stack 5 | 50 | ‘118 | | 1-221 -0150 


(600 stripped) 


Bristol 6 | 50 071 | -049| -052| -260| 
(400 u g.b.) 


wos os. 


| | | 
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experiment. It appears that results here are not significantly different from 
the rest of the measurements given in Fig.7 a, b. But the data in Table III 
show clearly that because spurious scattering in individual tracks increases 
with cell length (¢) almost like 7*'/*, it cannot be distinguished from true 
multiple coulomb scattering by comparing results for various cell lengths. 


Since the increase of spurious scattering with cell length does not permit 
us to distinguish true from spurious scattering, we need additional informa- 
tion about a particle before we can determine, or obtain limits on, the 
product of its velocity and momentum from scattering measurements. 
There are three commonly occurring cases of additional information. 


It may be known from range or variation of ionization density along the 
track that a particle was slow and that the expected true multiple coulomb 
scattering should be much higher than scattering values measured on known 
high energy tracks; we can then neglect spurious scattering and make the 
approximation 4... ~ Miue. This approximation applies in many experi- 
ments which have been carried out with nuclear emulsions in the past, in 


particular in the experiments designed to determine the mass of slow K- 
mesons. 


It may be found that there is a second track parallel and close to the 
track under investigation. We shall only discuss the simplest and most common 
case, namely when the two tracks are due to particles of equal mass and 
energy. This applies to tracks found in cyclotron exposed emulsions and 
it applies to the showers of relativistic a-particles which represent the evapora- 
tion of a partially destroyed heavy primary nucleus of great energy. Such 
a case was illustrated in Fig. 1. We have then the following relations: 


— D, = V/2 (€r,67 + (11) 
(D, — D,*) = V2 ER,G (7) 
therefore 
= 4 [(D, — D.)? — (D, — D,*)*] (12) 
D* as before refers to co-ordinate readings displaced by 20, with respect 


to those which determine the second difference D. The subscript | and 2 
refer to the two tracks, and 4; contains only true and no spurious scattering. 


The well-known result, namely that relative scattering on closely spaced 


tracks is more reliable than measurements on individual tracks, is examplified 
by the tracks whose contour was given in Fig. 1, The measured values arg 
listed in Table IV, 
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TABLE IV 
Individual and Relative Scattering Measurements on Two Parallel 
a-Particle Tracks 
(Expressed in microns) 


Cell length 250 500 ju 1000 2000 
D, 248 -540 -937 1-31 
V2 
A, -237 53 94 1-31 
D. +545 -994 1-33 
De -258 -562 1-02 1-32 
066 
Vv 2 
As 239 54 99 1-32 
A, 048 -080 0715 072 


Table IV shows that the scattering measurements are highly reproduce- 
able even after shifting the track by 20. It also shows that the two tracks 
give identical results, proving again that reading errors and grain noise 
are unimportant. 


The relative scattering 4; of the two tracks does not exceed the noise 
value for cell lengths as high as 2,000. The corresponding energy of the 
a-particles taking account of statistical uncertainties is > 20 BeV/nucleon 
in fair agreement with the energy of 15 BeV/nucleon estimated from the 
angular spread of the a-particle fragments emanating from this collision. 
On the other hand the scattering value, obtained on each track separately 
(which includes spurious scattering), corresponds to very low energies. It 
should be emphasized again that this event occurs in an emulsion which 
throughout has fairly moderate distortion (< 50 covans). 


The third case of additional a@ priori information applies to all tracks 
studied in our experiment; we know that because of our selection criteria, 
true scattering must be much smaller than the observed scattering values. 
Therefore : 


Anus. 
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We can, therefore, get some information on the nature of the displace- 
ments responsible for spurious scattering by considering the statistical 
distribution of observed D-values. We shall not in this paper make a 
detailed investigation of the type of displacement which is responsible for 
the effect of spurious scattering. Such an investigation is most conveniently 
carried out with plates which are exposed in the Bevatron and will contain 
many parallel tracks of mono-energetic high energy particles. Some tentative 
information on the nature of the dislocations has, however, been deduced 
from our data in the Appendix. 


TV. Spurious SCATTERING IN ITS RELATION TO THE CHARGE 
AND ENERGY SPECTRUM OF PRIMARY COSMIC RAYS 


We have shown that the determination of multiple coulomb scattering 
in emulsions on tracks of high energy particles is strongly affected by small- 
scale dislocations which occur during processing and whose effect on tracks 
is indistinguishable from true scattering. 


At present the range of validity of scattering measurements is limited 
by this spurious scattering and not, as has been generally believed, by 
various other types of noise and errors of measurement. Statistically as 
can be seen from Fig. 7, the spurious scattering parameters a’j),, determined 
for many tracks, have a broad distribution ranging from low values up to 
a maximum at a fairly sharply defined upper limit. This limit lies between 
-025 and -030°; it decreases slightly with increasing cell length and increases 
with increasing dip angles above ~2°. The occasional occurrence of 
very low scattering values does not contradict these findings. Spurious 
scattering values appreciably larger than the statistically determined upper 
limit do also occur occasionally. 


Spurious scattering is in our opinion responsible for the marked differ- 
ence between the energy spectra of primary nuclei obtained by Dainton, 
Fowler and Kent,® by Gottstein,’ and by Waddington!* all of whom used 
ordinary scattering measurements and the spectra obtained by Kaplon, 
Peters, Reynolds and Ritson® who used re/ative scattering measurements, 
nuclear interactions, the latitude effect and the range of heavy nuclei in 
absorbers. The measurements which we have discussed in the previous 
section show that only a small fraction of the relativistic primaries can be 
identified as relativistic by ordinary scattering measurements. 


For instance, there is general agreement about the total number of 
Primary nuclei with atomic number Z > 6 incident at high latitude (geo- 
magnetic latitude 54-55°). But there is no agreement on what fraction of 
these nuclei carry energy in excess of 1-6 BeV/nucleon. Experiments not 
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involving scattering measurements are consistent with each other; they 
include investigations with nuclear emulsions® !? rocket-borne ionization 
chambers,® cloud chambers triggered by a proportional counter telescope,’ 
and counter telescopes which include a Czerenkov detector.* These experi- 
ments yield a flux of heavy particles with energy « > 1-6 BeV/nucleon which 


is about three times as large as the flux determined from scattering measure- 
ments in emulsion.® 1% 


The charge determination on heavy primary nuclei can also not be carried 
out safely in experiments which involve multiple scattering measurements 
in emulsion. Table V shows the percentage of high energy particle tracks 
for which spurious scattering is so large (a’39, > *0125°) as to lead to an 
underestimate of charge by one or more units. The percentage of wrong 


charge assignments has been given separately for tracks of various inclina- 
tion and for various cell lengths. 


TABLE V 
Dip angle Percentage of underestimated charges derived 
cotg @ from scattering measurements 
Cell 250 500 750 1000 
lengths 

8 to 12 64% 43% 60% 15% 

12 to 20 50% 47% 39% 30% 

20 to 30 60% 42% 48% 40% 

>30 47% 22% 25% 25% 


In the experiment of Dainton, Fowler and Kent, who attempted to 
determine the primary charge distribution from measurements of scattering 
and of 5-ray density, nearly 80% of the tracks fall into dip angle interval 
12 < cotg 9 < 20 and were scattered with 250, to 500, cells. Other 
investigators have used still steeper tracks and shorter cell lengths for similar 
measurements. Figure 7a and Table V show that under these conditions 
reliable charge determinations cannot be obtained, and that the charge will 
be underestimated in 50-60% of the cases. 


In the presence of a strong Carbon peak in the primary radiation this 
“must lead to a large spurious peak of Boron nuclei, and thereby to an 
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over-estimate of the relative intensity of the light elements Li, Be, B compared 
to the C, N, O group. 


Apart from the presence of spurious scattering in the measurements, 
there are other grounds on which the conclusion of these authors, namely 
that the primary intensity of Li, Be, and B is comparable to that of C,N 
and O, can be criticised. The number of Lithium nuclei was corrected for 
an assumed lack of detection efficiency and was raised from the observed 
value of 45 to the corrected value of 130. The detection efficiency had not, 
however, been measured directly but was inferred from the observed Boron 
flux, which in turn, as we have seen, must to a large extent be spurious. 


Another reason for the large primary Li, Be, B flux obtained by these 
authors, lies in the method adopted in correcting for the Li, Be, B nuclei 
produced by fragmentation of heavier primaries in the upper atmosphere. 
The fragmentation probability was taken from measurements carried out 
with relativistic particles,! but it was applied to the mostly non-relativistic 
nuclei incident at high latitude. Since it is known that the mean free path 
for nucleons in nuclear matter varies by factors of order 2 or 3 from © 
energies of 0-3 BeV/nucleon to very high energies (the maximum occurring 
near -6 BeV),!® it seems reasonable to expect different fragmentation prob- 
abilities at different latitudes. By studying the fragmentation of heavy 
primaries in substances of low atomic weight it has been shown’ that the 
probability for obtaining Li, Be, B fragments is fairly high at medium and 
high latitudes and that the number of Li, Be, B nuclei produced by frag- 


mentation at balloon altitudes is much larger than was assumed by Dainton 
et al. 


Since the charge distribution cannot be reliably determined by scattering 
measurements, there remains at present only one method which has been 
suggested for studying the charge distribution in the low energy part of the 
spectrum. This method consists in accurate measurements of ionization 
versus range and is confined to tracks of particles which actually come to 
rest in nuclear emulsion; it has been applied by Waldeskog,”° and has so 
far given no evidence for a large Li, Be, B flux. But more work is needed 
to improve statistics in this experiment. 


Measurements in the high energy part of the spectrum are easier to 
make and methods are more numerous. 


Bradt and Peters, who used a stack of glass backed emulsions with 
alternating high and low sensitivity, obtained at geomagnetic latitude 
. A= 30° as an upper limit for the ratio Ry between the primary flux of Lithium, 
Beryllium, Boron nuclei and that of Carbon, Nitrogen and Oxygen nuclei 


= 
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at the top of the atmosphere the value Ry,<0-1. Their experiment has 
recently been repeated and extended to include primaries of charge Z = 2.! 
This last investigation includes a check of detection efficiency, which was 
found to exceed 85% for tracks with grain density corresponding to four 
times minimum ionization (primary a-particles) and must, therefore, be 
very close to 100% for tracks due to particles of charge Z> 3. 


Stix? in two very beautiful and difficult experiments identified particles 
at balloon height and geomagnetic latitude 4 = 41° by the pulse height in 
three proportional counters and at the same time by tracks in a cloud 
chamber triggered by these counters. He obtains Ry = 0-05+°3?. The 
upper limit is calculated by including dubious cases as belonging to the 
Li, Be, B group and by using not the standard deviation but the 95% prob- 
ability for the result to fall within the limits. 


Webber® at the same latitude, used a Czerenkov counter within the 
solid angle of a Geiger counter telescope. He obtains at balloon altitude 
a ratio LiBe B/Z>6= 0-16. Correction for light nuclei produced by 
fragmentation above the balloon will reduce this ratio, but the statistical 
errors are large and uncertain. 


Existing evidence at latitudes up to 41° favours, therefore, a low intensity 
of Li, Be and B nuclei. Since it is known that interstellar gas contains 
mostly hydrogen and that collisions of heavy nuclei (Z > 6) with hydrogen 
lead to disintegrations in which frequently the incident nucleus suffers only 
partial destruction, some Li, Be, B must be created in such collisions and 
must, therefore, be present in the primary beam. The observed intensity 
sets an upper limit to the fraction of heavy primaries which have suffered 
collisions in interstellar space and, therefore, to the amount of matter which 
they have traversed since their acceleration. 


In order to calculate this amount of matter one has to consider sepa- 
rately the fragmentation of the C, N, O group, and that of heavier primaries, 
because the processes by which Lithium, Beryllium and Boron are produced 
in collision with hydrogen are basically different, in the two cases.22 The 
evaporation of a nuclear fragment of charge Z> 3 in a collision between 
a heavy nucleus and hydrogen is rare, but its complete disintegration into 
singly and doubly charged particles is also very rare and the most common 
collision results in a reduction of charge of the heavy nucleus. On the 
other hand the collision of C, N, O nuclei with hydrogen must lead to the 
production of Li, Be and B residues except in those cases where fragmenta- 
tion into singly and doubly charged fragments is complete. 
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The probability (p) of producing a Li, Be, B fragment in a collision 
between a nucleon and a carbon, nitrogen or oxygen nucleus has been 
measured by various investigators. 


Harding®? by applying charge conservation to the prongs of cosmic 
ray stars in gelatine found that in 70% of the collisions the residue must have 
atomic number Z> 2. (Hodgson®* who carried out an experiment 
similar to Harding’s, obtained a value for p which is certainly much 
too low because it cannot be corrected for the large number of one and two 
pronged stars missed in his experiment.) 


The method of measuring the probability p used by Noon, Kaplon and 
Ritson seems preferable for our purpose because it is based on collisions 
of the same type as those presumed to occur in interstellar space. Out of 
87 nuclear interactions produced in an emulsion block detector by the heavy 
primary component incident at geomagnetic latitude A = 41°, they selected 
those due to collisions with protons characterised by the absence of evapora- 
tion tracks, by a recoil proton at much larger angle than the fast fragments 
of the incident nucleus, and by the fact that charge balance can be obtained 
on the assumption that the collision occurred with hydrogen in the target 
material. They obtained: : 


p=0-8+0°5 


It would be very desirable to reduce the statistical error by further observa- 
tions. Their result, however, is in rough agreement with that of Harding. 


If we assume p = 0-7 and remember that the group comprising C, N, O 
nuclei dominates the production by collision of Li, Be, B, one obtains an 
upper limit for the amount of interstellar hydrogen gas traversed by the 
primary cosmic radiation equal to 14% of the interstellar mean free path of 
C, N, O nuclei. We may take 300 millibarns as the cross-section for the 
interaction of fast proton with nitrogen and obtain then an upper limit of 
0-75 g./cm.? for the amount of interstellar material which the primary 
cosmic radiation has traversed in the time interval which elapsed between 
acceleration and arrival at the outer atmosphere. In view of the statistical 
uncertainty in the parameter p a more refined calculation does not as yet 
seem justified. This upper limit corresponds to the amount of matter 
encountered in approximately a thousand traversals of the thickness of the 
galactic disc. 


One can obtain a different upper limit by the following considerations. 
Iron is the heaviest element definitely identified among the primary particles. 
Iron nuclei constitute 4% of the primary radiation in a given velocity interval 
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(at a given geomagnetic latitude the percentage is larger) and carry 2% of 
the nucleons. If one makes the extreme assumption that at the time of 
acceleration the primary beam contained only iron nuclei one obtains from 
the fact that 2% did not suffer destructive collisions with interstellar hydrogen 
an upper limit of four mean free paths or 9-1 g./cm.? of hydrogen. If on 
the other hand one makes the more plausible assumption that at the time 
of acceleration the ratio of Iron to Carbon nuclei was the same as the 
abundance ratio of these elements in the universe (approximately 0-2) one 
obtains an upper limit of 0-7 mean free paths or 1-6 g./cm.? of hydrogen, 
not much higher than that derived from the abundance minimum in the 
primary Lithium, Beryllium and Boron component. 


Thus, it seems safe to conclude from the chemical composition of pri- 


mary nuclei, that few if any of the primary particles approaching our atmo- 
sphere have suffered collisions in the interstellar space. 
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APPENDIX I 


Figure 1 shows clearly that the magnitude of the dislocations in processed 
emulsions often varies appreciably in distances as small as 250u. If the 
direction of the displacements also changes within such small distance, one 
might expect to obtain a gaussian distribution G (D) for the individual values 
of the second differences D. On the other hand if the direction of displace- 
ments persist over larger distances comparable with the length of the tracks 
themselves, one must expect a quite different distribution. In that case, 
the resulting D-values will be proportional to the sine of the angle 4 between 
the track and the displacement vector and if this angle is arbitrary, the 
distribution of D-values takes the form: 


. 
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Here x, represents the contribution of ordinary noise to the spread of D- 
values and has been taken as x) = 0-075. In Fig. 8, we have plotted a _histo- 
gram of 826 D-values obtained with 1,000, cells on tracks whose dip angle 
is less than 1-9°. The gaussian distribution and the distribution F (D) have 
also been plotted. The distribution curves are normalized to the same area 
as the histogram and the parameter x,/x, = 12-3 has been chosen so as to 
yield the experimental value for the mean of the absolute second differ- 
ence, namely D = 0-513 p. 

Since F (D) fits the experimental distribution while G (D) does not fit, 
we may tentatively conclude, that the dislocations which give rise to spurious 
scattering change in magnitude over short distances but persist in direction 
for much longer distances. If this conclusion is correct, we should find in 
the same emulsion tracks with low and tracks with high spurious scattering. 


This seems to be borne out by our data. In Fig. 9, we have plotted 
the energies per nucleon obtained with 1000, cells for some tracks which 
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traverse several emulsion sheets. The graph indicates that energy values 
for the same particle are quite erratic from plate to plate, and that a single 
emulsion contains in general tracks which have large spurious scattering 
as well as tracks which have comparatively low spurious scattering. This 
we interpret tentatively as connected with the angle between the track and 
the dislocation vectors. 


Finally, if the spurious scattering is proportional to the sine of the angle 
% between track and dislocation vector, and if these angles are arbitrarily 
distributed, one expects for a large number of tracks a certain probability 
distribution for values of a’, namely: 


2 da 


where a) is determined by the observed upper limit of a’j9. We have only 
17 tracks with dip angle cotg @ > 30 which have been scattered with 1000, 
cell length (Fig. 75) and can, therefore, make only a rough comparison 


between expected and observed distribution of spurious scattering values. 
We find: 


Predicted Observed 
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